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ABSTRACT 

The X-ray emission from the wind-wind coUision in short-period massive 0+0-star 
binaries is investigated. The emission is calculated from three-dimensional hydrody- 
namical models which incorporate gravity, the driving of the winds, orbital motion of 
the stars, and radiative cooling of the shocked plasma. Changes in the amount of stellar 
occultation and circumstellar attenuation introduce phase-dependent X-ray variabil- 
ity in systems with circular orbits, while strong variations in the intrinsic emission 
also occur in systems with eccentric orbits. The X-ray emission in eccentric systems 
can display strong hysteresis, with the emission softer after periastron than at corre- 
sponding orbital phases prior to periastron, reflecting the physical state of the shocked 
plasma at these times. 

Our simulated X-ray lightcurves bear many similarities to observed lightcurves. 
In systems with circular orbits the lightcurves show two minima per orbit, which are 
identical (although not symmetric) if the winds are identical. The maxima in the 
lightcurves are produced near quadrature, with a phase delay introduced due to the 
aberration and curvature of the wind collision region. Circular systems with unequal 
winds produce minima of different depths and duration. In systems with eccentric 
orbits the maxima in the lightcurves may show a very sharp peak (depending on the 
orientation of the observer), followed by a precipitous drop due to absorption and/or 
cooling. We show that the rise to maximum does not necessarily follow a 1/dsep law. 
Our models further demonstrate that the effective circumstellar column can be highly 
energy dependent. Therefore, spectral fits which assume energy independent column(s) 
are overly simplified and may compromise the interpretation of observed data. 

To better understand observational analyzes of such systems we apply Chandra 
and Suzaku response files, plus poisson noise, to the spectra calculated from our simu- 
lations and fit these using standard XSPEC models. We find that the recovered tem- 
peratures from two or three-temperature mekal fits are comparable to those from fits 
to the emission from real systems with similar stellar and orbital parameters/nature. 
We also find that when the global abundance is thawed in the spectral fits, sub-solar 
values are exclusively returned, despite the calculations using solar values as input. 
This highlights the problem of fitting oversimplified models to data, and of course is 
of wider significance than just the work presented here. 

Further insight into the nature of the stellar winds and the WCR in particular 
systems will require dedicated hydrodynamical modelling, the results of which will 
follow in due course. 

Key v^fords: shock waves - stars: binaries: general - stars: early-type - stars: mass 
loss - stars: winds, outfiows - X-rays: stars 



1 INTRODUCTION 

The shock-heated plasma in the wind-wind collision 
region (WCR) of massive stellar binaries can produce 
copious X-ray emission. The emission often displays 
* E-mail: jmp@ast.leeds.ac.uk orbital variability, which can result from changes to 
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the occultation of the emitting region by the stars, 

to the attenuation through the stellar winds, and to 

I i I 

t he separat io n of the stars (e.g. iRauw e t al. 2002 



Schild et al, 

20061: iLinder et alj 



2004 _ ISana et alJ |2004|: iDe Beck e r ct al] 
20061 : IPollock fc CorcoranI l2006l: 



Naze et all l2007l : iHamaeuchi et all l2007l : ISana et al] 
20081 ). In wider systems, the post-shock plasma may 



exhib it signs of non-equilibrium ionization (|Pollock et alJ 
l2005l). and non-equilibri um electron and ion temperatures 
IZhekov fc Skinneij[2OO0l ). 

The X-ray emission is a useful probe of the underlying 
wind parameters. The hardness of the emission is related to 
the post-shock temperatures within the WCR, which in turn 
depends on the pre-shock wind speed. The X-ray bright- 
ness depends on the pre-shock density of the winds, while 
the absorption of soft X-rays through the circumstellar en- 
vironment depends on t he integrated density along sight 
lines to the WCR fe.g. IStevens. Bloridin fc Pollock! 1 19921: 
Stevens et alj Il996l: fpittard fc Steyend Il997l: IPittard et alj 



19981: iPittard fc CorcoranI |2002|: iParkin fc PittardI l2008l: 



Parkin et ai]|2009l ). Hence both the brightness and the de- 



gree of absorption provide information about the stellar 
mass-loss rates. Although the X-ray emissivity is propor- 
tional to the square of the density, inhomogeneties can be 
rapidly smoothed out within adiabatic WCRs: thus the re- 
sulting X-ray emissi on may be re latively insensitive to the 
presence of clumps (|Pittardll200"7l ). Since mass-loss rate es- 
timates are often uncertain due to unknown wind clumping 
factors, an insensitivity to clumping potentially allows the 
X-ray emission from CWBs to provide accurate determina- 
tions of stellar mass-loss rates. Recent observations in the 
UV have highlighted the uncertainty which still exists, with 
mass-loss rate estimates differi ng by factors of up to 100 with 
respect to other methods (e.g. Bouret. Lanz fc Hillier 20051 : 
iMartins et al] l2005bl : iFuUerton. Massa fc Prinial 120061 '). A 
recent review of the curre nt situation can be found in 
IPuls. Vink fc Naiarrd (|2008l l. 

The WCR can also be a site of particle acceleration. 
The energetic particles produce no n-thermal radio emis- 
sion via the synchr otron process (e.g. iDoughertv et al.ll2003l : 
IPittard et ani2006l ). and non-thermal X-ray and 7-ray emis- 
sion from inverse Compton c ooling, neutral pion decay , 
and relativistic bre msstrahlung (|Pittard fc Doughertvl[200g : 
iLevder et all l2008h . The non-thermal radio emission can 



DouEhertv. WiUiams fc PoUaccd l2000bl; 


IDoughertv et alJ 


20051). and can also undergo dramatic variations in flux (e.g. 


WiUiams et al.l 19921: White fc Becker 


1995: 


Rauw et alJ 


2002b: De Becker et al.1 l2004d: iBlomme et alj 


20051. I2OO7I: 



Van Loo et al.ll2008l ). If the particle acceleration efficiency 
is sufficiently high, the thermal structure of the WCR may 
be affected, resulting in softer X-ray emission as the plasma 
becomes cooler and denser. In this way the characteristics of 
the thermal X-ray emission may also constrain the efficiency 
of particle acceleration. 

Models of the X-ray emission from colliding wind sys- 
tems based on hydrodynamical simulations have, to date, 
been almost entirely performed in two-dimensions, with 
an underlying assumption of axissymmetry. While this ap- 
proach is perfectly reasonable for wide systems with long or- 
bital periods, axissymmetry is a poor assumption in shorter 
period systems where orbital effects become important. 



Though three-di mensiona l hydro dyn amical simulations hav e 
been presented bv lWalderl (119981 ) and lLemaster et all (|2007h . 
these works also assumed that the winds were instanta- 
neously accelerated to their terminal velocities. In reality, 
the winds in short-period systems collide prior to reaching 
their terminal velocities, so realistic simulations must also 
account for the acceleration of each wind. In a new advance, 
three-dimensional mo dels with radiat ively driven winds were 
recently presented by iPittardI (|2009al . hereafter Paper I). In 
addition to the acceleration of the winds, these models also 
account for orbital motion of the stars, gravity, and cooling 
in the post-shock plasma. 

In this work we examine the thermal X-ray properties 
of the WCR from these models. We produce synthetic X-ray 
spectra and lightcurves, and examine how these change with 
the viewing angle. Section [2] describes details of the models 
and summarizes the method of calculating the X-ray emis- 
sion and absorption. This section also notes the procedure 
adopted for folding the theoretical spectra through the re- 
sponse files of current X-ray observatories to simulate "fake" 
observations, which are subsequently fit using standard anal- 
ysis techniques. We present our results in Section |3l Com- 
parisons to previous numerical models and observations are 
made in Sections U and [Sj respectively. Section|B]summarizes 
and concludes this work. 



2 DETAILS OF THE CALCULATIONS 
2.1 The Numerical Models 

The X-ray calculations in this paper are based on the 
three-dimensional hydrodynamical models described in Pa- 
per I. The models incorporate t he radiative driving of 
the st ellar winds (based on the ICastor. Abbott fc Klein! 
(1 19751) formalism, with the finite disk correction factor 



Pauldrach. Puis fc Kudritzki! (|l986l )). gravity, orbital ef- 



of 

fects, and cooling. The models were not designed to be of 
particular systems. Rather, the aim was to obtain a bet- 
ter understanding of how the nature of the collision region 
depends on various key parameters. The models are summa- 
rized in Tables [Tj and [2] The assumption of main sequence 
stars minimizes the effects of tidal distortions, which are 
not modelled. The winds are also assumed to be spherically 
symmetric. Further details about the models can be found 
in Paper I. 

In model cwbl two identical 06V stars move around 
each other in a circular orbit with a period of 3 days. The 
stellar separation is 34.26 R0, and each star has an or- 
bital velocity Worb = 290 kms~^ . The thermal behaviour of 
the WCR can be described by the ratio of the cooling time 
to the characteristic flow time of the hot shocked plasma, 

X ~ ^, where vg, is the pre-shock wind speed in units 

of 1000 km s^^ , d\2 is the separation of the stars, and Af_7 
is the stellar mass- loss rate in units of lO~^M0yr~^ (c.f. 
IStevens et al]|l992! ). In model cwbl, the WCR is highly 
radiative (x << 1), and significantly distorted by orbital 
effects, showing strong aberration and downstream curva- 
ture. Model cwbl is simil ar to many real sys tems, including 
HD 215835 (DHCep; seelLinder et al.1 1200?!. and reference s 



therein), HD 16 5052 (lArias et al 
and HD 159176 jPe Becker et al ' 



l2002l: 



2004b! : 



Linder et al] | 2007|) . 



Linder et aL!l2007l ) 
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Table 1. Assumed binary parameters for the models calculated in Paper I. The semi-major axis is 34.26 R0 in model cwbl, 76.3 Rq 
in models cwb2 and cwb3, and 55 R© in model cwb4. e is the orbital eccentricity and r) is the (terminal velocity) momentum ratio of 
the winds, forb and ''w are the orbital speeds of the stars and the preshock wind speeds along the line of centres, x is the ratio of the 
cooling time to the characteristic flow time of the hot shocked plasma. X ^ 1 indicates that the shocked gas rapidly cools, while X 1 
indicates that the plasma in the WCR remains hot as it flows out of the system. Larger values of the ratio iiorb/^w produce a greater 
aberration angle, 9^^^, and tighter downstream curvature, of the WCR. The degree of downstream curvature of the WCR in the orbital 
plane is given by Qcoriolisi where the curvature is assumed to trace an Archimedean spiral which in polar coordinates is described by 
^ = c«coriolisS- The value of Ocoriolis corresponds to the approximate downstream distance (in units of ciscp) along the WCR for each 
radian of arc it sweeps out in the orbital plane. Smaller values indicate tighter curvature. The leading and trailing arms of the WCR 
in model cwbS display differing degrees of curvature, so the value quoted for this model is an average. The pre-shock orbital and wind 
speeds in model cwb3 are also different for each star /wind - the flrst (second) value is for the primary (secondary) star/wind. The values 
of Xi ^orb/'^w, Ssib S'lid Cfcoriolis phase dependent in model cwb4, because of its eccentric orbit - values at periastron and apastron 
are quoted. The values for Cfcoriolis a^rc calculated after comparing the orbital speeds at periastron and apastron against those in models 
cwbl and cwb2, and represent the "instantaneous" curvature at these phases. 



Model 


Stars 


Period 


e 


V 


forb 




X 


^orb/fw 


^ab 


*-^coriolis 






(d) 






( km s"-"- ) 


( km s^-^ ) 






n 


(dscp rad^-"-) 


cwbl 


06V+06V 


3 


0.0 


1 


290 


730 


0.34 


0.40 


17 


3.5 


cwb2 


06V+06V 


10 


0.0 


1 


225 


1630 


19 


0.14 


3-4 


6.5 


cwb3 


06V+08V 


10.74 


0.0 


0.4 


152,208 


1800,1270 


28,14 


0.16 - 0.084 


~ 2 


4.5 


cwb4 


06V+06V 


6.1 


0.36 


1 


334 - 156 


710 - 1665 


0.34 - 19 


0.47 - 0.09 


21-4 


3-10 



All of these systems have near identical main-sequence stars 
of spectral type 06— 07, and circular or near-circular orbits 
with periods near 3 days. 

In model cwb2 the orbital period is increased to 10 
days, with the stellar separation becoming 76.3 Rq. The 
winds collide at significantly higher speeds than in model 
cwbl, and the postshock gas is largely adiabatic. The aber- 
ration and downstream curvature of the WCR are both less- 
ened compared to model cwbl. Model cwb2 is similar to 
HD 93161 A, an 08V + 09V system with a circular orbit 
and an orbital period of 8.566 days (|Naze et ahlboOSl ). al- 
beit with slightly more massive stars and powerful winds. 
Another system which is not too dissimila r in its properties 
is Plaskett's star (|Linder et al.ll2006l . [ioosh . though this ob- 
ject contains stars which have evolved off the main sequence. 

Model cwb3 examines the interaction of unequal winds 
in a hypothetical 06V-I-08V binary. The stars in this model 
have the same separation as those in model cwb2. The pri- 
mary wind collides at higher speed than the secondary wind, 
and its postshock plasma is slightly more adiabatic. 

Model cwb4 investigates the effect of an eccentric orbit 
which takes the stars through a separation of 34.26—76.3 Rq 
(i.e. the separations of the stars in the circular orbits of mod- 
els cwbl and cwb2). The WCR is radiative at periastron and 
adiabatic at apastron, and the aberration and downstream 
curvature are phase dependent. A surprising finding from 
Paper I is that dense cold clumps formed in the WCR at 
periastron persist near the apex of the WCR until almost 
apastron. This is because the clumps have relatively high 
inertia, and fiow out of the system much more slowly than 
the hotter gas which streams past them. Some well-known 
O-l-O binaries with eccentric orbits include (in order of in- 
creas ing orbital period) HP 152248 (e = 0-127: ISana et al] 
I2OO4I ). HD 9 3205 (e = 46: iMorreU et al.ll2001f ). HD 93403 
(e = 0.234; iRauw et al .'"2002a'), C ygOB2#8A (e = 0.24 ; 
iDe Becker. Rauw fc Ma nfroid 2 0oi: iDe Becker et all 120061 ) 
and tOrionis (e = 0.764: iBagnuolo et al.ll200ll ). 



2.2 Modelling the X-ray emission and absorption 

To calculate the X-ray emission we read our hydrodynam- 
ical models into a radiative transfer ray-tracing code, and 
calculate appropriate emission and absorption coefficients 
for each cell using the temperature and density values. A 
synthetic image on the plane of the sky is then generated 
by solving the radiative transfer equation along suitable 
lines of sight through the grid. Since non-equilibrium ef- 
fects are small in short period O-l-O systems (see Paper I), 
the X -ray emissivity is c alculated using the mekal emission 
code (|Mewe et al.l 1 19951 . and references therein) for opti- 
cally thin thermal plasm a in collisional ionization equilib- 
rium. Solar abundances (| Anders fc Grevessd Il989l ) are as- 
summed throughout this work. The emissivity is stored in 
look-up tables containing 200 logarithmic energy bins be- 
tween 0.1 — 10 keV, and 91 logarithmic temperature bins 
between 10* - 10^ K. Line emission dominates the cooling 
at temperatures below 10^ K, with thermal bremsstrahlung 
dominating at higher temperatures. The hydrodynamical 
grid is large enough to capture the majority of the X-ray 
emission from each of the models. 

The main contributors to the absorption of keV X-rays 
are the K shells of the CNO element s. The photoe lectric 
absorption is calculated using Cloudy (|Ferlandll20dol ). The 
opacity is stored in look-up tables containing 26 temper- 
atures between 10* — lO'' K. As in previous works (e.g. 



Luo. McCrav fc Mac Lowlll990l:IStevens. Blondin fc PoUockl 



19921 : iMvasnikov fc Zhekovlll993l : iPittard fc Stevend [l997l ) . 
electron scattering is neglected. Electron scattering becomes 
important once the optical depth to this process nears unity, 
i.e. when To — A^eO"e ^ 1, where A^o is the column density of 
free electrons along a line of sight and is the Thomson 
cross-section. In the ionized winds, the proton and electron 
column densities are approximately equal. In our models, 
Tc < 1 is indeed satisified. For example, in the dense cir- 
cumstellar environment of model cwbl, an observer viewing 
the system pole on (i = 0°) sees an average "effective" hy- 
drogen column density to high energy (2 — 10 keV) X-rays of 
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2 X 10^^ cm"2 (see Fig.[5l;a) and Section [3T3l for further 
details). Since occulation is minimal, the "effective" column 
density in this case reflects the true column density of the 
circumstellar environment. The electron scattering optical 
depth is then Tc ~ 0.13. Note that the higher "effective" 
column densities shown in Fig. [3a) for an observer in the 
orbital plane {i = 90°) at phase 0.0 are instead a reflec- 
tion of the occultation that takes place at this time, and do 
not indicate that electron scattering becomes optically thick 
(see Section 13.1.31 for further details). Electron scattering 
will, however, be important in systems with denser winds, 
where very high column densities can be reached. Such sys- 
tems, inc lude, for example, the supermassive system rj Car 
(see, e.g.. |Parkin et al.ll2009l ). The likely effect is that abrupt 
changes in the emission (e.g. in lightcurves and spectra) will 
be somewhat smoothed/blurred, though we leave a study of 
this effect to future work. 

The present calculations also have an interstellar ab- 
sorption column of 10^^ cm"^ added to them, and each 
model is assumed to be at a distance of 1 kpc from an ob- 
server. The X-ray spectra/lightcurves were calculated from a 
single "frame" (i.e. changing the orientation) for the circular 
orbit models (cwbl, cwb2, and cwb3), and from a sequence 
of snapshots for the eccentric model (cwb4). The lightcurve 
for i = 0° is invariant for models cwbl, cwb2, and cwb3. 



Table 2. Assumed stellar parameters for the models. 



Parameter / Star 


06V 


08V 


Mass (Mq) 


30 


22 


Radius ( ) 


10 


8.5 


Effective temperature (K) 


38000 


34000 


Mass-loss rate (Mq yr^^ ) 


2 X 10-7 


10-7 


Terminal wind speed (kms^^ ) 


2500 


2000 



The spectra were binned with the FTOOLS task grp- 
pha so that each energy bin contained a minimum of 20 
counts. The "fake" spectra were fitted using XSPEC version 
12.5.0ac, distributed with HEASoft6.6.3. Since an ISM col- 
umn of 10^^ cm-^ was added to our theoretical spectra, we 
force the absorbing column to each model component to be 
at least as large. However, we note that if this restriction is 
relaxed, the best fit often returned lower columns. The the- 
oretical spectra were generated using emissivities from the 
mekal thermal emission code, so for consistency we also fit 
the data in XSPEC using the mekal thermal model. 



2.3 Generating and analyzing "fake" spectra 

In the following section we "observe" the theoretical spec- 
tra generated from our ray-tracing code with the Chan- 
dra and Suzaku X-ray observatories. This involves convolv- 
ing the theoretical spectra with the energy response and 
effective area of these telescopes, to generate spectra in 
counts/energy bin. Counting statistics are included in this 
process. The resulting "fake" spectra are then analyzed us- 
ing XSPEC, and fitted with standard spectral models, in 
an analogous manner to the analysis of real data (the only 
difference is that a background component does not need to 
be subtracted). The aim is to study how the fit parameters 
compare with those from the analysis of real data, and how 
they compare to what is known about the theoretical input 
spectra. This type of analysis remains very novel, having 
been applied to colliding wind binarie s only by Stevens ct al.l 
(fl996), Pittard fc Stcvci- i s (199 7), Z hekov fc Skin ner (2003), 
and lPittard fc Corcoran! (|2002l ). 

The majority of our analysis is concentrated on sim- 
ulated Suzaku XIS spectra. To generate these we used the 
XISO ancillary response file (ARF) and redistribution matrix 
file (RMF) for an on-axis point source downloaded from the 
HEASARC websitsQ- Although these are old (2006) calibra- 
tions, they are fine for our purpose, which is to investigate 
the values and variation of the best-fit parameters, and the 
corresponding fiuxes of the best-fit models. A small number 
of simulated Chandra ACIS-I spectra were also computed. 
These used the Cycle 11 ACIS-I aimpoint ARF and RMF, 
downloaded from the Chandra websit^EI. 



^ http:/ /heasarc. gsfc.nasa.gov/docs/heasarc/caldb/data/ 
suzaku /xis/ index, html 

^ http:/ /cxc. harvard.edu/caldb/prop_plan/imaging/indcx. html 



3 RESULTS 
3.1 Model cwbl 

3.1.1 Images 

Fig.[T]shows broad-band images from model cwbl for an ob- 
server located directly above the orbital plane (i = 0°). The 
stars are oriented north-south in these images (cf. the images 
of the thermal radio emission in |Pittard.2009bl , hereafter Pa- 
per II) , and the orbital induced aberration and downstream 
curvature of the WCR is clearly visible. The emission mor- 
phology reflects the underlying structure and dumpiness of 
the WCR, resulting from the powerful dynamical instabili- 
ties present in this system. A detailed discussion of the hy- 
drodynamics can be found in Paper I. The projected emis- 
sion from different inhomogeneities merges together near the 
apex of the WCR, but individual clumps and bowshocks can 
be identified further downstream. It is clear that a small 
amount of emission, particularly at the lowest energies, is 
not captured due to the finite extent of the hydrodynamical 
grid used in the model, but this loss should not be signifi- 
cant. The hard (2 — lOkeV) emission predominantly arises 
from the apex of the WCR. Although there are regions of 
hot gas further downstream (see Paper I) , the density there 
is too low for these regions to contribute significantly to the 
emission. The spatial scale of the emission is far too small to 
be resolved with current X-ray telescopes: WR 147 is likely 
to remain the only s ystem with a spatially resolved WCR 
(IPittard et al. I l2002d ) for some time to come. 

Fig.[2]shows broad-band images from model cwbl for an 
observer located in the orbital plane {i = 90°). The clumpy 
nature of the WCR and the bowshocks around some of the 
denser regions are visible. At particular phases/ viewing an- 
gles the emission from bright parts of the WCR is occulted 
by the foreground star. Additional foreground emission is 
sometimes seen in front of the stellar disc at these moments. 
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Figure 1. Broad-band intensity images from model cwbl at i = 0°. From left to right the images are in the 0.1 — 0.5 keV, 0.5 — 2.5 keV, 
and 2.5 — lOkeV bands. The grayscale covers 4 orders of magnitude in brightness, with black corresponding to a maximum intensity 
of 10^ ergcm"^ keV~^ ster~^ in the left and right panels, and 10* ergcm"^ keV~^ ster~^ in the central panel. The major ticks 
on each axis mark out 0.2 mas. The stars are at ±0.08 mas north and south of the image centre. The arrows mark the direction in 
which an observer looks into the system at the indicated orbital phases. We also define an azimuthal viewing angle, <j>, which increases 
anti-clockwise from the bottom of the image, so that at phase 0.0 4> = 0° . </) = 90° corresponds to phase 0.75, and = 180° and </) = 270° 
correspond to phases 0.5 and 0.25, respectively. 
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Figure 2. Broad-band 0.5-2.5 keV intensity images from model cwbl at j = 90°. From left to right the phase of the observation increases 
from 0.0 (conjunction, (j> = 0°), to 0.125 (</> = 315°), to 0.25 (quadrature, (p = 270°), to 0.375 {<f> = 225°). The grayscale covers 4 orders 
of magnitude in brightness, with black corresponding to a maximum intensity of 10* erg cm" -2 s-l keV-l ster-l. The major ticks on each 
axis mark out 0.2 mas. 
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Figure 3. X-ray lightcurves for models cwbl (left), cwb2 (middle), and cwb3 (right), calculated over the energy bands 0.1 — 0.5 kcV 
(top), 0.5 — 2.5 keV (middle), and 2.5 — 10 keV (bottom) for inclination angles i = 30°, 60° and 90°. In all cases, the observer is located 
along a direction vector specified hy (j) = 0° (i.e. the longitude of pcriastron, to = 270°). The stars are at conjunction at phases 0.0 and 
0.5, and quadrature at phases 0.25 and 0.75. 
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3.1.2 Lightcurves 

Lightcurves from model cwbl, computed over the energy 
bands 0.1 - 0.5 keV, 0.5 - 2.5 keV, and 2.5 - 10 keV, 
are shown in the left column of Fig. [S] The stars pass 
in front of each other at phases 0.0 and 0.5, and are at 
quadrature at phases 0.25 and 0.75. As the orbit is cir- 
cular, the intrinsic emission is constant, so the variations 
displayed in the lightcurves in Fig. [3] are entirely due to 
changes in the occultation and circumstellar absorption as 
a function of phase. If there were no orbital induced ef- 
fects on the WCR (i.e. no aberration or curvature of the 
WCR), the lightcurves would display dual symmetry about 
phases corresponding to both quadrature (0.25, 0.75) and 
conjunction (0, 0.5) of the stars to the line of sight (c. f. 
iPittard fc Stevens|[l997l : lAntokhin. Qwocki fc Brownll2004 ) , 
because of the equal winds and constant stellar separa- 
tion. When orbital effects are included, the symmetry about 
quadrature is broken, and the lightcurves instead are ex- 
pected to show a double periodicity in the case of identical 
winds. However, close examination of Fig. [3] reveals that 
in fact this dual periodicty is also broken (note that the 
peaks in the 0.1 — 0.5 keV lightcurve have different heights). 
Clearly, the dynamical instabilities which form in the WCR 
develop independently in each arm, and break this symme- 
try too. 

The variation of flux with orbital phase is largest in all 
lightcurves when the inclination angle i = 90° , and decreases 
with decreasing i. The amplitude of variation is also largest 
in the softest band. Both of these findings are expected: the 
soft emission is more easily absorbed by the stellar winds, 
while lines of sight to the WCR pass, on average, through a 
greater attenuation column, and there is also greater occul- 
tation of the emission region by the stars, when the observer 
is in the orbital plane. 

Minima in the lightcurves occur near phases 0.0 and 
0.5, when the emission from the WCR suffers the great- 
est reduction by stellar occultation and wind absorption. 
However, close examination reveals that the minima actu- 
ally occur slightly before each conjunction. This refiects the 
aberration of the WCR. In model cwbl the abberation angle 
is ~ 17°, which corresponds to 0.05 in orbital phase, and is 
therefore similar to the observed lead. 

Short, sharp dips are also seen in the lightcurves near 
orbital phases 0.32 and 0.82 (being most visible in the 
0.5 - 2.5 keV lightcurve for i = 90°). This is due to ab- 
sorption from the thin dense layer of cooled post-shock gas 
(s ee Paper I). The dips are similar to those seen in Fig. 15 
of lAntokhin et all !\2004 ). but are broader and less obvious 
due to the orbital-induced curvature of the WCR. They also 
display a phase lead consistent with the phase lead of the 
main minima. 

The ISM corrected 0.5 - lOkeV luminosity is 6.75 x 
10^^ ergs"^ 1.41 x 10^^ ergs~^, and 6.54 x 10^^ ergs"^ 
at viewing angles of {i,4>) = (90°, 0°), (90°, 90°), and (0,0) 
(pole-on), giving Lx/iboi = 4.7 x 10"^, 9.9 x 10"^, and 4.6 x 
10~^, respectively. These values are all significantly above 
th e scaling law (log L x/Ijboi)= —6.912 ± 0.153) determined 
by ISana et all (|2004l ). and are thus indicative of a binary 
system with strong colliding winds emission. 



3.1.3 Spectra 

The emission from model cwbl is overall very soft, reflecting 
the rapid cooling of the hot plasma in the WCR. Fig. |3Ia) 
displays synthetic X-ray spectra as a function of orbital 
phase for an observer in the orbital plane. The spectra at 
phase 0.5 and 0.75 are almost identical to those at phase 0.0 
and 0.25, differing only due to the time dependence of the 
dynamic instabilities present in each arm of the WCR, and 
so are not shown. The hard emission is lower at conjunction 
(phases 0.0 and 0.5) when one of the stars passes in front 
of the apex of the WCR, than at quadrature (phases 0.25 
and 0.75). This is a combination of occultation and wind ab- 
sorption, with the former dominating at hard energies (see 
below). The softer emission is also reduced at these phases 
as it is attenuated by a greater amount of unshocked stel- 
lar wind between the WCR and the observer. In contrast, at 
quadrature the emission suffers less attenuation and occulta- 
tion, and it more easily escapes the system. Fig.lJJb) shows 
that the variation of the soft emission with the inclination 
angle, i, is small. The soft emission arises from relatively low 
temperature plasma, comprising a relatively large volume of 
the WCR. Consequently the loss of flux due to circumstellar 
absorption and stellar occultation is relatively independent 
of the orientation. As we shall see, interstellar absorption is 
a large contribution to the overall absorption of soft X-rays. 

The variation of the hard emission with i (see Fig. |4}3) 
is of a similar magnitude to the change between conjunction 
and quadrature when the observer is in the orbital plane 
(see Fig. |4^) . This is not surprising since in both cases the 
variation is mainly due to changes in the occultation, and 
these changes are comparable: more of the apex of the WCR 
is revealed as the observer's sight line moves away from the 
eclipse at i = 90° , (p = 0° . 

Figs. [5ja) and (b) show the "effective" circumstellar 
column and optical depth as a function of energy for an ob- 
server in the orbital plane. These were calculated by compar- 
ing the ray-traced intrinsic and attenuated spectra (the lat- 
ter prior to the addition of the interstellar column). Their ra- 
tio gives an energy dependent optical depth (Fig. [5]3) , which 
can be converted into an "effective" column (Fig. [5^) by con- 
sidering the energy dependent opacity of the cold plasma in 
the winds. Note that this "effective" column does not always 
provide a good measure of the circumstellar column for an 
observed X-ray, since it can be substantially weighted by 
occultation effects (see below). Instead it is intended as an 
indicator of how realistic the usual method of using energy- 
independent columns with, for example mekal fits, is. 

In Fig. [5ja) we see that the circumstellar column from 
the phase 0.0 calculation monotonically increases with en- 
ergy (except at a few noteable absorption edges), reflecting 
the fact that the hardest emission is on average generated 
closest to the stars (see Fig. [1]). The effective attenuation 
resulting from circumstellar absorption and stellar occulta- 
tion exceeds that due to interstellar absorption (assumed to 
be Nn = 10^^ cm~^) at energies above 0.2 keV. The effec- 
tive circumstellar column exceeds 10^^ cm~'^ for E > 4keV. 
Columns this high cannot be produced by the winds them- 
selves, and instead indicate substantial occultation of the 
intrinsic emission by the stars. The strong energy depen- 
dence of the column at phase 0.0 as shown in Fig. [5ja) also 
hints at weaknesses that are inherent in simple analyses 
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Figure 4. (Top) X-ray spectra from models cwbl (left), cwb2 (middle), and cwb3 (right) for an observer in the orbital plane (i = 90°) 
as a function of orbital phase. (Bottom) As top but as a function of i for an observer with = 0°. 
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Figure 5. "Effective" circumstellar column (top) and optical depth (bottom) from models cwbl (left), cwb2 (middle), and cwb3 (right) 
for an observer in the orbital plane (i = 90°) as a function of energy and orbital phase. The circumstellar column and optical depth for 
an observer at j = 0° are also shown. Note that in some cases the "effective" column can be strongly weighted by occultation effects - 
at these times it is not a good indicator of the actual circumstellar column that the observed X-rays travel through. 



of X-ray data of colliding wind binaries in which energy- 
independent columns are applied in model fits. While in- 
dividual columns to multiple emission components clearly 
provide some level of energy-dependent A^h, this is of course 
achieved in a rather crude and clumsy way. 

The equivalent optical depth from circumstellar absorp- 
tion and occultation is shown in Fig. \5][h). It is highest at 
the lowest energies. At phase 0.0, the optical depth r = 6 
at 0.1 keV, dropping to r = 0.26 at lOkeV. In comparison. 



the ISM column provides an optical depth of 30 at 0.1 keV, 
declining to 10~^ at lOkeV. 

Occultation of the harder X-ray emission is much re- 
duced at phase 0.25 (see the spectrum in Fig. S^), and this 
is manifest as a marked drop in the "effective" circumstel- 
lar column and optical depth values in Figs. [5ja) and (b). 
In fact, the column is now roughly constant over the energy 
range 0.5 — 7 keV. In this case a simple analysis using energy- 
independent columns would probably not be too bad an ap- 
proximation at this phase. However, it is clear that such 
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analyses have major short-comings in short period binaries 
at phases when occultation is hkely to be significant (for 
instance, the intrinsic luminosity can be significantly under- 
estimated - see the following section). At the very highest en- 
ergies {E > 7keV) the effective circumstellar column again 
declines. This is due to the fact that the highest tempera- 
ture plasma exists in localized bowshocks around clumps far 
downstream in the WCR, as discussed in Paper I. 

Figs. I2a) and (b) also show the effective circum- 
stellar column and optical depth for an observer directly 
above/below the centre of mass in the orbital plane (i.e. at 
an inclination angle i — 0°). The degree of absorption and 
occultation are intermediate between the conjunction and 
quadrature phases of an observer in the orbital plane. 

The deduced values of the energy-dependent circumstel- 
lar column shown in Fig. [5Ta) can be co mpared to analytical 
estimates. Using Eq. 11 in lStevens et al. (1992,), which gives 
the column density at quadrature from the stagnation point 
(assumed to be on the line of centres in the axisymmet- 
ric case) through the undisturbed terminal speed wind, we 
obtain A'^h ~ 3.4 x 10^^ cm~^. While this is in good agree- 
ment with the column to the higher energy emission (which 
should be the best proxy to the apex of the WCR) at phase 
0.25 (see Fig. [S^), it would see m to be a som e what fortu- 
itous coincidence. For instance, IStevens et al.l (|l992l ) note 
that under the assumption of axisymmetry the column at 
quadrature is independent of the system inclination. How- 
ever, Fig. [5ja) shows that the circumstellar column for a 
pole-on observer (i — 0°) is significantly greater at most 
energies than the phase 0.25 column. This is because with 
equal strength winds all the material along the line-of-sight 
from the apex of the WCR to an observer at i = 0° has been 
processed through the WCR, and is denser than the sur- 
rounding unshocked winds. In contrast, the line-of-sight to 
an observer in the orbital plane at phase 0.25 passes mostly 
through unshocked wind material. So there are actually large 
differences in the density structure along th ese two sight- 
lines, yet this is not accounted for in Eq. 11 of lStevens et al.l 
(|1992D . 

3.1.4 Spectral fits 

In this section we "observe" the theoretical spectra gener- 
ated from model cwbl with the Chandra and Suzaku X-ray 
observatories. Table [3] notes the results of the subsequent 
spectral fits, where an exposure time of 10 ksec has been as- 
sumed. At least three mekal components are needed to ob- 
tain satisfactory fits to the simulated Suzaku spectra, while 
the simulated Chandra spectra require at least two mekal 
components. 

Fig. [nia) shows the fake Chandra spectrum from model 
cwbl at i = 90° and phase 0.0, with Poisson statistics added. 
A single temperature wabs(mekal) model is a very poor fit to 
the simulated data (xl = 3.91). However, an acceptable fit 
is achieved with the addition of another mekal component 
(xt = 1.27). Fig. |6jb) shows the corresponding spectrum 
and fit at phase 0.25. The hotter temperature component 
has a slightly reduced temperature, a higher absorbing col- 
umn, and a lower normalization at phase 0.0 (conjunction) 
compared to phase 0.25 (quadrature), in line with expecta- 
tions (see Figs. and [5^). 

Two mekal components return reasonable fits to sim- 



ulated Suzaku spectra with an exposure time of 10 ksec 
(xt = 1.37, at phase 0.0 when i = 90°), but fail to provide 
a good fit when the exposure time is increased to 40 ksec 
(xt ~ 1-64), notably failing to fit the line at 0.55 keV. 
Adding a third mekal component does not significantly im- 
prove the fit (xt = 1.50). Irrespective of the orbital phase, 
the temperatures of the three mekal components are all be- 
low 0.75 keV, reflecting the relatively cool shocked plasma 
created by the relatively low preshock wind speeds in this 
model. Furthermore, the best-fit absorbing columns are of- 
ten significantly higher than the ISM value, indicating that 
the circumstellar absorption of X-rays is significant in this 
model. 

The circumstellar columns returned from the spectral 
fits are consistent with the trend shown in Fig.[5ja) of higher 
effective column with energy. Also of note is that the derived 
A'h to the hottest mekal component of the spectral mod- 
els is greater at phase 0.0 (conjunction) than at phase 0.25 
(quadrature). In addition, the hottest mekal component is 
both hotter and brighter at phase 0.25. These are consistent 
with the changes to the circumstellar absorption and stellar 
occultation of the WCR with the orientation of the observer, 
as noted earlier. 

A good fit to the simulated Chandra spectrum can also 
be obtained with a single mekal component if the global 
metal abundance is allowed to vary. In this case, a metal 
abundance of 0.131'^q q2 ^<3 provides a good fit {xi = 1.11) 
to the data. This is very interesting, since we know that the 
actual plasma has solar abundances. Clearly, there is great 
opportunity for the analysis of low spectral resolution data 
to return unphysical fit parameters when fitting emission 
from inherently multi-temperature plasma with simpler (e.g. 
single temperature) spectral models. In an identical (10 ksec) 
analysis to a fake Suzaku spectrum we again find that a sin- 
gle temperature fit is very poor [xt = 2.69), but this time 
it remains poor [xt = 2.13) when the global metal abun- 
dance is allowed to vary (fitting at z — 0.079 ± 0.011 Zq). 
This nicely illustrates the advantage of having higher spec- 
tral resolution. The story is more complicated for fits with 
two mekal components. Tying the global abundances of the 
mekal components together, one finds that the returned 
value from the analysis of the phase 0.0 fake Chandra spec- 
trum is very poorly constrained, with its value depending on 
how the model approaches its minimum (e.g. the initial val- 
ues entered into the model). The 90 per cent confidence in- 
terval typically extends from metallicities of 0.4 — 2. Ox solar, 
with a "best-fit" value oi z = 0.79 Zq. A three-temperature 
fit to the fake Suzaku phase 0.25 spectrum returns a global 
abundance of z = 0.71 Zq (since xt > 2, the uncertainty 
on this value cannot be estimated using the "error" com- 
mand in XSPEC). It would therefore appear that the fits re- 
turn more accurate abundances the more complex the spec- 
tral model is. These finding s agree with the earlier work of 
[Strickland fc Stevend (|l998i ) who were examining ROSAT 
data in a related context. 

Table [4] shows the intrinsic luminosity of model cwbl, 
plus the intrinsic luminosity returned from the spectral fits 
(i.e. the observed luminosity, corrected for the interstellar 
and circumstellar absorption determined by the fit). We find 
that the intrinsic luminosities returned from the fits in all 
cases underestimate the true intrinsic luminosity from the 
system, by factors of up to 2. This discrepancy arises because 
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'fake" spectra from the models for different viewing angles, a) Two-temperature fit to a fake Chandra spectrum 
= 90° and phase 0.0 (conjunction), b) As a) except the system is viewed at quadrature (phase 0.25). c) Three- 
fake Chandra spectrum of model cwb2 for an observer in the orbital plane viewing the system at conjunction (phase 
to a fake Suzaku spectrum. Note the much better spectral resolution in the Suzaku spectrum. The spectra in panels 
"exposures" of lOksec. e) Three-temperature fit to a fake Suzaku spectrum of model cwb3 for an observer in the 
; the system at phase 0.0 (weaker 08V wind in front), f) As e) but at phase 0.5 (stronger 06V wind in front). Note 
strength of the Fe K line. The spectra in panels e) and f ) have effective "exposures" of 20 ksec. 



direct occultation of the emission is a significant factor in 
close binaries like cwbl, yet no account is made for occulta- 
tion in the fits. It again highlights problems which can ensue 
when fitting too-simple models to data. The discrepancy is 
larger at conjunction than at quadrature, as expected. We 
further note that the size of the discrepancy is not depen- 
dent on the resolution of the two spectra (i.e. Chandra ver- 
sus Suzaku). The discrepancy reduces in wider systems (see 
model cwb2, next), but could be even larger in yet closer 
systems. 

3.2 Model cwb2 

3.2.1 Images 

Fig. [7] shows the broad-band intensity images from model 
cwb2 for an observer located directly above/below the or- 
bital plane. The hardest emission is again confined to a re- 
gion close to the apex of the WCR, and the curved shape 
of the WCR is clearly seen. The leading edge of each arm 
of the WCR is sharper and more distinct. The trailing edge 
is blurred because the hot plasma inside the WCR displays 
an increasing phase lag as one goes further from the orbital 
plane. This reflects the fact that the motion of this plasma 
is strongly influenced by the prior (rather than the current) 
positions of the stars. The extent of the low surface bright- 
ness emission from the WCR is affected by the size of the 
numerical grid. Calculations with a bigger grid would re- 
veal that, for instance, the left-hand edge of the emission in 
the bottom right corner of the 0.5 — 2.5 keV image in Fig. [7] 
would extend further to the left. However, we remain con- 



fident that the majority of the emission is captured in this 
(and the other) models, since the surface brightness is 4 or- 
ders of magnitude lower than the peak surface brightness 
obtained at the apex of the WCR. 

Broad-band intensity images from model cwb2 for an 
observer in the orbital plane are shown in Fig. [8] The images 
bear substantial similarities to the radio images shown in 
Paper II. The morphology of these images is determined by 
the relative orientation of the observer to the WCR and the 
stars. The foreground star eclipses the emission from those 
parts of the WCR which lie behind it. The double-helix-like 
structure seen in the soft and medium-band images when the 
system is half-way between conjunction and quadrature (at 
phase 0.375) is due to limb-brightened emission. The vertical 
curvature again illustrates the increasing phase-lag of the 
shocked gas with distance above/below the orbital plane. 
The increasing confinement of emission to the apex of the 
WCR at higher energies means that the hard-band images 
instead show a disc-like structure. In all 3 energy bands the 
surface brightness of the emission is highest in the images 
generated at quadrature, when the X-rays from the apex 
of the WCR initially escape through the hot, low opacity, 
gas within the WCR. This is also true of radio images when 
u < 100 GHz (see Paper II). 

Images for an observer viewing the system with an incli- 
nation angle i — 30° are shown in Fig. |§1 The brightest parts 
of the WCR are again those which are limb-brightened. The 
overall morphology of the WCR is "S" -shaped, and the po- 
sition of the foreground star is again clear through its occul- 
tation of background emission. Comparison to radio images 
again reveals significant similarities (see Fig. 8 in Paper II). 



10 J. M. Pittard and E. R. Parkin 



10.1-0.5 keV 



1 1 0.5-2.5 keV 





12.5-10 keV 



Figure 7. Broad-band intensity images from model cwb2 at i = 0°. From left to right the images are in the 0.1 — 0.5 keV, 0.5 — 2.5 keV, 
and 2.5 — lOkeV bands. The maximum intensity in the images is 10^ ergcm"'^ s~^ keV~^ ster~^, and the grayscale covers 4 orders of 
magnitude in brightness. The major ticks on each axis mark out 0.5 mas. The stars are located north and south of the image centre by 
0.18 mas. 
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Figure 8. Broad-band intensity images from model cwb2 at i = 90°. From top to bottom the images are in the 0.1 — 0.5 keV, 0.5 — 2.5 keV, 
and 2.5 — lOkeV bands. From left to right the phase of the observation increase from 0.0 (conjunction, (p = 0°), to 0.125 (</> = 315°), to 
0.25 (quadrature, = 270°), to 0.375 {4> = 225°). The maximum intensity in the 0.1 - 0.5keV images is 10^ ergcm'^ s"'^ kcV"^ ster'^. 



while in the 0.5 — 2.5 keV, and 2.5 — lOkeV images it is 10^ erg cm 
brightness, and the major ticks on each axis mark out 0.5 mas. 
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Figure 9. Broad-band 0.5 — 2.5keV intensity images from model cwb2 at i = 30°. From left to right the orbital phases are 
0.0, 0.125, 0.25, and 0.375. The stars and the WCR again rotate anti-clockwise, while the maximum intensity of the images is 
10^ erg cm~^ s~^ keV~^ ster~^ , and the major ticks on each axis mark out 0.5 mas. 
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3.2.2 Lightcurves 

Lightcurves from model cwb2 are shown in the central col- 
umn of Fig.|3l The longer orbital period and the higher speed 
of the wind-wind collision changes the lightcurves compared 
to those from model cwbl in several important ways: 

(i) The luminosity in the 2.5 — lOkeV lightcurve is over 
an order of magnitude higher. This reflects the harder emis- 
sion resulting from the higher postshock temperatures cre- 
ated by the faster wind collision speeds (1630 kms~^ versus 
730 kms"^ at the apex of the WCR). 

(ii) The amplitude of variation with orbital phase is much 
reduced in the 0.5 — 2.5 keV and 2.5 — lOkeV lightcurves. 
This reflects weaker circumstellar attenuation due to the 
lower wind densities surrounding the WCR, and reduced 
occultation effects due to the larger size of the WCR relative 
to the stars. 

(iii) The lightcurves display clear symmetry, with 2 iden- 
tical cycles per orbit. Dynamical instabilities are weak in 
this model, because the WCR is largely adiabatic and the 
winds have equal speeds, and do not appreciably disrupt the 
inherent symmetry between each arm of the WCR. 

(iv) The sharp absorption features seen in the lightcurves 
from model cwbl have disappeared, since there is no longer 
a dense thin layer of post-shock gas to absorb the X-rays in 
this way. 

The ISM corrected 0.5 — lOkeV luminosity is 1.55 x 
10^^ ergs-^ 1.73 x 10^^ ergs'^ and 1.77 x 10^=* ergs"^ 
at viewing angles of (i, 0) = (90°, 0°), (90°, 90°), and (0,0) 
(pole-on), giving Lx/Lboi = 1.09 x 10"^ 1.21 x 10~^ and 
1.24 x 10~^, respectively. The luminosities and Lx/iboi val- 
ues are slightly greater than from model cwbl, and are again 
consistent with strong colliding winds emission. 

3.2.3 Spectra 

X-ray spectra from model cwb2 are shown in Figs. |4jc) 
and (d). It has already been noted that the continuum emis- 
sion from model cwb2 is much harder than from model cwbl, 
due to the higher pre-shock speeds which the winds attain 
before they collide in this system. However, the line emis- 
sion also reflects the higher temperatures in model cwb2: 
a strong FeK line is visible at approximately 6.7 keV, but 
this line is much weaker relative to the continuum in model 
cwbl. In contrast, most of the other spectral lines in model 
cwb2, particularly those with E < 1 keV, display weaker 
emission relative to the continuum than in model cwbl, re- 
flecting the lack of strong cooling in model cwb2 compared 
to model cwbl. There is significantly more soft emission at 
1 = 0° than at higher inclination, which reflects the fact 
that the sight line to the stagnation point of the WCR is 
entirely through the hot lower opacity WCR. This enhance- 
ment does not occur in model cwbl, due to the different 
nature of the WCR (specifically the relative lack of hot gas 
within it) in this model. The soft and hard emission is again 
lower at conjunction due respectively to enhanced circum- 
stellar absorption and occultation, as was also the case for 
model cwbl. 

Figs. Efc) and (d) show the "effective" circumstellar 
column and optical depth as a function of energy from 
the ray-traced calculation for a variety of viewing angles. 



The column at conjunction (phase 0.0) shows again a gen- 
eral rise with increasing energy. The columns (and hence 
the optical depths) at 0.1 and lOkeV are similar to those 
from model cwbl, though they are substantially reduced in 
value at intermediate energies (e.g. at 1 keV the circumstel- 
lar column and optical depth is now an order of magni- 
tude lower). Lower circumstellar columns and optical depths 
are expected, of course, because of the wider stellar separa- 
tion. In model cwb2 the circumstellar column at phase 0.25 
peaks at around 1 keV, whereas in model cwbl the peak 
column at phase 0.25 is to emission near 4keV (ignoring 
the high column to line emission at 6.7 keV). Fig. ^[c) also 
shows that the circumstellar column to an observer with 
i = 0° is lower than the columns obtained for an observer 
in the orbital plane, consistent with the higher flux of soft 
X-rays at this orientation (see Fig. |4ji) . This contrasts with 
model cwbl where the column for a pole-on system lies (for 
E > 0.25 keV) inbetween the columns to observers with 
i = 90° at phase 0.0 and 0.25. 

Another noticeable difference to model cwbl is that the 
effective columns show significant variability between adja- 
cent energy bins. This is caused by the different tempera- 
tures (and thus locations) at which line emission and the ad- 
jacent continuum are formed. Finally, we can again make a 
compariso n between the c ircum stellar column obtained from 
Eq. 11 in IStevens et all (|l992l ') and Fig. ^c). The former 
gives A^H ~ 1.5 X 10^^ cm~'^, which is substantially greater 
than the values shown in Fig. [5jc) to an observer in the or- 
bital plane at quadrature (phase 0.25) and to an observer 
with i — Q° . This further highlights that the formula in 
I Stevens et al] (|l992l ') is unsuitable for use in short-period 
systems. We therefore make no further comparisons to it in 
this work. 

3.2.4 Spectral fits 

Fig. inic) shows a simulated 10 ksec Chandra spectrum from 
model cwb2 at i = 90° and phase 0.0. Folding the same theo- 
retical spectrum through the Suzaku response and exposing 
for 10 ksec yields the "fake" spectrum shown in Fig. [S^d) . 
The higher spectral resolution of the latter observatory is 
clearly evident. Two-temperature fits to these spectra are 
uniformally poor (typically xt ^ 2), with the flux at low 
and high energies underestimated. Another mekal compo- 
nent is clearly required. As expected three-temperature fits 
are more acceptable. Importantly the fits to both the Chan- 
dra and Suzaku spectra return consistent parameter values. 
At phase 0.0, both sets of fits find that the normalization 
of the mekal components increases monotonically with the 
temperature of the component (see Table [Sj. Compared to 
the results from model cwbl, the returned temperatures are 
significantly higher, and the absorbing columns significantly 
lower, both of which are consistent with expectations. 

Another finding is that the mekal components consis- 
tently return cooler temperatures at phase 0.25 compared 
to phase 0.0. This reflects the greater ease at which low 
energy X-rays can escape absorption by the circumstellar 
environment at this phase (see Fig.|4l:), and is further man- 
ifest by the consistently lower columns to the mekal com- 
ponents at phase 0.25 compared to phase 0.0. In all four 
of the flts (to the "fake" Chandra and Suzaku spectra at 
phases 0.0 and 0.25) the absorbing column to the lowest 
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Figure 10. Broad-band 0.5 — 2.5 keV intensity images from model 
cwb3. The left panel shows the emission at i = 0°. The right panel 
displays the emission at i = 30° and <j) = 45° (phase 0.875). 
The maximum intensity is 10^ ergcm"^ s~^ keV~^ ster~^, and 
the major ticks on each axis mark out 0.5 mas. ij> = 0° corre- 
sponds to the 08V star in front. 

temperature mekal component (at kT ~ 0.32 keV) is in- 
dicative of only ISM absorption. Fig. [5jc) reveals that the 
additional circumstellar column is indeed small in com- 
parison (« 2 X 10'^° cm~^, versus 10^^ cm~^ for the as- 
sumed ISM column). The additional (above ISM) columns 
to the medium and hard mekal components returned from 
the fit to the phase 0.0 Chandra spectrum (respectively 
1.5l'i':5 X 10^^ cm"^ and 2.ltl l x 10^^ cm'^), though not 
particularly well constrained, are comparable to the columns 
shown in Fig. [5jc) obtained at the energy of the individual 
mekal components. This is a pleasing result. Surprisingly, 
although the implied circumstellar columns to the medium 
mekal components obtained from the phase 0.0 and 0.25 
Suzaku fits are also comparable to those shown in Fig. [5jc), 
the fits find that there is no need for additional circumstellar 
absorption to the hard mekal components. In this respect, 
the fits to the "fake" Chandra spectra are better at recover- 
ing the actual circumstellar columns than the higher spectral 
resolution Suzaku spectra. Having said this, we note that the 
columns to the hard component from the Suzaku fits have 
upper limits which are still consistent with Fig. [SJc) . 

We again note that when fitting a three-temperature 
mekal model to medium-resolution spectra, relaxing the 
global abundance in the model can lead to erroneous results. 
For the i — 90° phase 0.0 case, the fit to the Chandra spec- 
trum returns z = O.GOj^JJ jg Zq. From the Suzaku spectrum 
we obtain z — 0.73^^']^ ^Q- Both results are significantly 
below the solar abundances used in our models. Finally, we 
again compare the real intrinsic luminosities calculated di- 
rectly from the models and the inferred intrinsic luminosities 
from the spectral fits. Table |4] shows that while luminosity 
differences still exist, occultation effects are now largely in- 
significant. Indeed, the returned intrinsic luminosity is now 
often greater than the intrinsic value. Such differences, in- 
cluding overestimates, result from the imperfect nature of 
the fit and poisson noise in the count rate. 

3.3 Model cwb3 

3. 3. 1 Images 

The stellar winds in model cwb3 are of unequal strength, 
being blown from hypothetical 06V and 08V stars. The 
stronger wind from the 06V star pushes the WCR closer 
towards the 08V star, and, compared to model cwb2 where 



the winds are of equal strength, bends the arms of the WCR 
inwards towards the weaker wind. The left panel of Fig. 1101 
shows broad-band 0.5 — 2.5 keV intensity images from model 
cwbS for an observer directly above/below the orbital plane. 
A comparison against the corresponding image in Fig. [7] re- 
veals several differences. Firstly, the brightest part of the 
image (at the apex of the WCR) is located closer to the 
08V star (which is to the south in these images). Secondly, 
the downstream positions of the arms of the WCR are also in 
different locations, due to the different ram-pressure balance 
of the winds. Thirdly, there are differences in the brightness 
contrast across the contact discontinuity. At the apex of the 
WCR, ram pressure balance requirements mean that since 
the 08V wind has a lower pre-shock velocity than the 06V 
wind, the 08V material must have a higher pre-shock den- 
sity. This directly translates into a higher post-shock density 
(and lower post-shock temperature), and thus into greater 
X-ray emissivity on the 08V side of the contact discontinu- 
ity. A similar effect is also seen from the thermal radio emis- 
sion (see Paper II). This emission contrast is then amplified 
in the leading arm of the WCR, but reduced in the trailing 
arm, as a result of the dynamics of the WCR (see Paper I). 
In comparison, there is initially no contrast in the emission 
across the contact discontinuity at the apex of the WCR in 
model cwb2, though such an effect subsequently develops in 
the downstream arms. Finally, we note that due to velocity 
shear Kelvin-Helmholtz instabilities occur along the contact 
discontinuity in model cwbS, and these are also visible in 
the images. The clearest sign occurs at phases 0.375 and 
0.875 in Fig. 1111 where an oscillation of the contact discon- 
tinuity separating the bright and fainter parts of the WCR 
can be seen. The right panel in Fig. [7] shows the intensity 
image for an observer with i = 30° and = 45°. The 08V 
star is silhouetted against the WCR. The emission from the 
shocked 08V wind is clearly brighter, for the reasons already 
described above. 

Fig. [Tl]shows broad-band 0.5 — 2.5 keV intensity images 
from model cwb3 for an observer in the orbital plane with 
= 0° as a function of phase. The 08V star is in front 
at phase 0.0, while the larger 06V star is in front at phase 
0.5. Significant differences to the images from model cwb2 
are again apparent (cf. Fig. [S| . For instance, at conjunction 
when the weaker 08V wind is in front (phase 0.0), the limb 
brightened edge of the leading arm of the WCR (to the right 
side of the image) shows greater curvature and is projected 
closer to the centre of the image, while the limb brightened 
part of the trailing arm is located at the far left side of the 
image. In addition, the 08V star occults a smaller region of 
the WCR than the larger 06V star does in model cwb2. 

Other differences are also apparent. At phase 0.375, the 
double-helix-like structure seen from model cwb2 is replaced 
by a more imbalanced morphology, where the brightest re- 
gions trace the hmb-brightened edge of the dense, shocked 
08V gas on the trailing edge of the leading arm. Interest- 
ingly, compared to simulated radio images (see Fig. 12 in 
Paper II), at phase 0.875 the shocked 06V gas is much more 
visible on the right side of the image. Obviously the size of 
the occulted region differs depending on which star is in 
front, but we note that the position of this relative to the 
limb-brightened part of the WCR is also different at phase 
0.375 and 0.875, reflecting the different amounts of down- 
stream curvature imparted to the leading and trailing arms 
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Figure 11. Broad-band 0.5 — 2.5 keV intensity images from model cwb3 at i = 90°. From top to bottom and left to right the orbital 
phase is 0.0 (0 = 0°, conjunction, 08V star in front), 0.125, 0.25 = 270°, quadrature), 0.375, 0.5 {<j> = 180°, conjunction, 06V star in 



front), 0.625, 0.75 ( 
each axis mark out 0.5 mas 



90°, quadrature), and 0.875. The maximum intensity is lO^crgcm ^s ^ kcV 



- ster 



and the major ticks on 



of the WCR (see Paper I for further details). Finally, we note 
that the vertical curvature of the WCR is also apparent in 
the images at quadrature (phase 0.25 and 0.75). 



3.3.2 Lightcurves 

The right column of Fig. O shows the orbital phase variation 
of the X-ray emission from model cwb3. The 08V star is 
in front at phase 0.0, with the 06V star in front at phase 
0.5. The unequal wind strengths in model cwb3 is manifest 
in the unequal depths of the two minima around the orbit. 
The deeper minimum occurs when the 06V star (which has 
the denser wind) is in front of the WCR apex, as expected. 
A phase lead to the bottom of the minimum is present in 
some lightcurves (e.g. the 0.5 — 2.5 keV curve), though not in 
others (e.g. the 2.5 — 10 keV curve) . The emisison at quadra- 
ture (phases 0.25 and 0.75) is brighter than at conjunction, 
as seen in the other simulations. The luminosity in the soft- 
est lightcurve is slightly higher near phase 0.25 than near 
0.75. This is because lines of sight to the apex at phases 
near 0.25 initially pass through the hot, low opacity, gas in 
the WCR because of orbital aberration. 

The X-ray luminosity is somewhat lower in this model 
compared to model cwb2, reflecting the reduced wind power 
of the 08V star. The reduction in luminosity is greatest 
in the harder energy bands (the 2.5 — lOkeV luminosity 
in model cwb3 is only 40 per cent of the luminosity in 
model cwb2), reflecting the slower speed of the 08V wind 
and the increased obliquity of the shock in the primary 
wind with off-axis distance relative to model cwb2. The 
ISM corrected 0.5 - 10 keV luminosity is 1.07 x 10^^ erg s~^ , 
1.14 X 10^^ ergs~^ , and 1.16 x 10^^ ergs~^ at viewing an- 
gles of [i, (j)) = (90°, 0°), (90°, 90°), and (0, 0) (pole-on). The 
slightly reduced X-ray luminosities compared to model cwb2 
are somewhat offset by a corresponding reduction in Lboi, 
giving //x/iboi ~ I0~^ for all 3 of the viewing angles indi- 
cated, again consistent with strong colliding winds emission. 
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Figure 12. Intrinsic X-ray spectra from each shocked wind in 
model cwbS. 



3.3.3 Spectra 

The right column of Fig. |4] shows intrinsic and attenuated 
spectra from model cwb3. Because of the unequal winds, the 
emission at phases 0.5 (0.75) is no longer identical to that at 
phases 0.0 (0.25), so all conjunction and quadrature phases 
are shown. There is now a clear difference in the strength 
of the low energy absorption at orbital phases 0.0 and 0.5, 
reflecting changes to the wind density along sight lines when 
either the 08V or 06V stars are in front. 

In models cwbl and cwb2, the winds were of equal 
strength and the intrinsic emission from their shocked 
plasma was identical, and contributed equally to the to- 
tal. However, Fig. [12] shows that the intrinsic emission from 
the postshock wind of the 06V star is harder than the in- 
trinsic emission from the postshock wind of the 08V star. 
This reflects the higher velocity at which the 06V wind en- 
counters the WCR, due to the greater distance over which 
it can accelerate and its higher terminal velocity (although 
this increase in the preshock velocity is somewhat offset by 
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the greater shock obhquity further downstream). But while 
the shocked 06V wind dominates the hard emission, the 
shocked 08V wind dominates the overall intrinsic emission, 
as is clearly apparent in Fig. 1121 contributing approximately 
two-thirds to the total luminosity. This is despite the total 
wind power of the 08V star being just one third of that of 
the 06V star. The fact that it dominates the X-ray luminos- 
ity is due to a greater fraction of it being processed through 
the WCR, and it subsequently rad iating energy more effi- 
ciently (c.f. lPittard fc Stevensli2002l 'l. 

Figs. EJe) and (f ) show the "effective" circumstellar col- 
umn and optical depth as a function of energy for various 
viewing angles into model cwb3. We find similar columns at 
phases 0.0 and 0.5 to those reported previously from models 
cwbl and cwb2. The column at phase 0.5 (when the larger 
06V star and its denser wind is in front) is slightly higher 
than at phase 0.0 (when the smaller 08V star is in front) 
as expected. The slope of the column with energy and its 
bin-to-bin variation are similar to those from model cwb2 
shown in Fig. [5jc). As previously noted, the different wind 
strengths in model cwb3 break the symmetry that models 
cwbl and cwb2 display at quadrature. Fig. [5je) shows that 
the columns at phase 0.75 are very low at high energies, 
whereas they are significantly higher at phase 0.25. In con- 
trast, the columns at iJ < 0.25 keV at phase 0.0 and 0.75 
are identical, whereas they are lower at phase 0.25. 

The effective circumstellar column to an observer at 
i = 0° is 1.7 X 10^^ cm~^ at 0.1 keV, rises an order of mag- 
nitude to ~ 3 X 10'^'^ cm~'^ by 1 keV, and then drops below 
10^^ cm~^ by 7.5 keV. The i = 0° circumstellar column and 
optical depth are greater than those from model cwb2, since 
the unequal wind strengths bend the WCR around the 08V 
star, so that only a relatively short segment along the line-of- 
sight from the apex of the WCR passes through hot plasma, 
with the rest of it through wind material from the 06V star. 

3.3.4 Spectral fits 

For model cwb3 we simulate Suzaku spectra with an expo- 
sure time of 20ksec. Two-temperature mekal fits are poor 
(xi ~ 2.0), and underestimate the hard X-ray flux. Three- 
temperature fits are much more acceptable (see Fig.[S^e) and 
(f), and Table|31). Significantly higher absorption is found to 
the hot component at phase 0.5 compared to phase 0.0, con- 
sistent with expectations given that the denser 06V wind 
is in front at phase 0.5. The temperature returned to the 
hottest mekal component is also significantly lower at phase 
0.5 compared to phase 0.0. This reflects the much weaker 
FeK line emission at phase 0.5, due to greater occultation 
(by a larger star) of the hottest part of the WCR at phase 
0.5 (the surface brightness of the high energy emission falls 
off very rapidly from the apex of the WCR, and thus the ob- 
served flux of high energy X-rays is highly susceptible to the 
amount of occultation). No extra absorbing column (above 
the ISM value) is required to the separate mekal components 
in many of the fits. This is roughly consistent with Fig. [5jc) 
where it can be seen that the effective circumstellar column 
is typically low, due to the relatively small mass-loss rates 
and wide stellar separations in the model. 

Having said this, significant additional absorption above 
the ISM value is returned from the fit made to the phase 0.5 
spectrum. The degree of extra absorption increases with the 



temperature of the mekal component (no extra absorption is 
required to the lowest temperature mekal component in the 
fit). The additional absorption to the 0.85 keV component 
(4tf X lO^^cm-^) 

is a bit lower than the effective circum- 
stellar column at this energy in Fig. [5je) (though within the 
90 per cent confidence range the fitted value is reasonable) . 
In contrast, the additional absorption to the 1.66 keV com- 
ponent (3.6^}'5 X 10^^ cm~^) is both in better agreement 
with Fig. \Sle) and is more tightly constrained. 

Interestingly, if the exposure time is reduced to 10 ksec, 
the fit to the i — 90° phase 0.0 spectrum returns fcTa — 
1.68lQ^g keV. Although formally this is still consistent with 
the value returned from the fit to the 20 ksec exposure spec- 
trum (fcTa = 2.09to 4o keV), it highlights that a lack of 
counts at high energies due to relatively short exposures 
can bias the resulting fits towards lower temperatures. 

3.4 Model cwb4 

3.4.1 Images 

In contrast to the previous models which all had circular or- 
bits, model cwb4 simulates a CWB with an eccentric orbit 
(e — 0.36). This introduces a time-dependence to the intrin- 
sic emission, which now varies with phase, whereas it was 
constant in the circular orbit models cwbl— cwb3. Fig. 1131 
shows intensity images of the 0.5 — 2.5 keV broad-band emis- 
sion from model cwb4 for an observer directly above the 
orbital plane. The images show striking variations in their 
brightness and morphology as a function of orbital phase, re- 
flecting the dramatic changes in the WCR during the orbit 
(see Paper I for full details of the hydrodynamics). Dramatic 
changes were also seen in the thermal radio-to-sub-mm emis- 
sion (see Paper II for further details). 

There is a smooth morphology to the emission from the 
WCR at phase 0.5 (apastron), when the WCR is adiabatic 
and instabilities are rare. As the stars progress in phase the 
WCR rotates in the images. The WCR becomes brighter and 
shows increasing curvature as the snapshots move towards 
the time of periastron passage (the maximum surface bright- 
ness of 1.9 X 10* ergcm^'^ s^^ keV~^ ster~^ occurs at phase 
0.95. At periastron there is a distinct change in the morphol- 
ogy of the images, with instabilities now clearly visible. This 
reflects the sudden cooling and formation of dense clumps 
within the WCR (see Paper I for further details). This mor- 
phology persists until phase 0.2, at which point the ratio of 
the cooling time to the flow time of the shocked gas near 
the apex of the WCR becomes signiflcant again. Dense cold 
clumps which formed during the periastron passage are still 
present in the WCR, but there is now also a substantial vol- 
ume of hot gas. These cold clumps are gradually destroyed 
or cleared out of the system, so that by phase 0.7 there is 
no longer any cold, post-shock, gas on the hydrodynamical 
grid (note, however, that even at phase 0.5 this process is far 
from complete) . Careful examination reveals that individual 
clumps and their ablated tails are visible in these images, al- 
though the level of detail is not great enough for this to be 
apparent in Fig. [131 Since the WCR rotates to follow the 
motion of the stars, while the dense clumps flow out on al- 
most ballistic trajectories, the clumps are often seen exiting 
the WCR through its trailing shock. They are then exposed 
to the full "fury" of whichever high-speed wind they find 
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Figure 14. As Fig. ll3l but showing broad-band (2.5 — 10 licV) intensity images. The grey scale and tick marks are the same as in Fig. 1131 
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themselves in, and are enveloped by a bowshock and high 
temperature plasma. In Paper I it was speculated that this 
additional contribution to the overall amount of hot plasma 
in the system may be significant in terms of the resulting 
X-ray luminosity. However, it is now clear from Fig. 1 131 that 
such regions have a negligible effect in this regard (although 
careful inspection reveals that there are in fact noticeable 
features in the intensity images). 

Fig. [14] shows intensity images of the 2.5 — 10 keV 
emission. The emission is clearly less extended and is 
more concentrated towards the central part of the WCR. 
Otherwise the morphology and brightness of the emis- 
sion behaves in a rather similar way to that in the 0.5 — 
2.5 keV images. The highest surface brightness of 3.5 x 
10^ erg cm~^ s"'^ keV~^ ster^^ occurs at phase 0.9. 

3.4-2 Lightcurves and spectra 

Fig. [15] displays the lightcurves which are obtained from 
model cwb4. The top panels show lightcurves of the intrinsic 
emission (red curves), while the bottom panels show atten- 
uated lightcurves. The intrinsic X-ray emission is no longer 
constant with phase, but, at low energies, reaches a max- 
imum at or near periastron as the stars reach their point 
of closest approach (Fig. I15b 'l. In contrast the 2.5 — lOkeV 
intrinsic emission (Fig. \15k ) peaks at phase 0.9, and then 
undergoes a precipitous drop so that it is actually close to 
its minimum value at periastron. 

The phase-dependence of the intrinsic emission is com- 
plicated, and depends on a combination of factors, including 
the current separation of the stars, their separation in the 
recent past, and the variation of the pre-shock wind densities 
and speeds and the post-shock cooling efficiency over this in- 
terval. The intrinsic emission from systems w here the WCR 
is lar gely adiabatic should scale as 1/dsop ([Stevens et al.l 
1 19921 ). This relationship breaks if cooling within the WCR 
becomes important (cf. Fig. 9 in [Pittard fc Stevens] [l997l ) . 
and /or if the pre-shock speeds of the winds change. Both 
of these events occur in model cwb4, so therefore it is not 
surprising that the intrinsic luminosity does not follow this 
scaling. In fact the intrinsic luminosity varies by factors 
of approximately 20, 4 and 20 in the 0.1 — 0.5, 0.5 — 2.5 
and 2.5 — lOkeV lightcurves, respectively, whereas for com- 
parison a 1/dsep response predicts only a factor of 2 vari- 
ation. For an observer with viewing angles of i = 90° 
and = 0°, the ISM corrected 0.5 — lOkeV luminos- 
ity varies between 1.22 x 10"^'^ ergs"^ at phase 0.1, and 
4.53 X 10'^'^ ergs~^ at phase 0.95, representing a change in 
ix/iboi from 8.56 x 10"^ to 3.18 x 10~®, respectively. For an 
observer located pole-on (i, </> = 0°,0°), Lx/iboi varies be- 
tween 7.01 X 10"^ at phase 0.1 to 5.88 x 10~® at phase 0.95. 
The peak values of I/x/iboi are the highest obtained from 
any of the models in this work, and indicate not only the 
strength of the colliding winds emission in this system, but 
also the relative ease with which the X-ray photons escape 
the system at favourable orientations because of the nature 
of the WCR. 

To understand the behaviour of the intrinsic lightcurves 
of model cwb4, one needs to examine the intrinsic spectra, 
which display strong phase-locked variation. Fig. 1161 shows 
that the spectral hardness of the intrinsic X-ray emission 
changes by a huge amount between periastron and apas- 



tron. At periastron the emission is very soft, reflecting the 
low pro-shock wind speeds at this phase (710kms~^ along 
the line of centres), whereas the emission is much harder at 
apastron, since the increase in the stellar separation allows 
the winds to accelerate to higher speeds before their colli- 
sion (1665 kms~^ along the line of centres), and reduces the 
effects of radiative inhibition (see Paper I). 

Interestingly, the intrinsic X-ray spectrum is at its hard- 
est at phase 0.6, rather than at apastron. Thereafter, the 
intrinsic emission begins to soften as the stars move closer 
together, the preshock wind speeds decline, and the WCR 
becomes increasingly radiative (see Fig. I17b l . The softening 
is initially manifest as an increase in the soft emission, while 
the harder emission (e.g. E > 4 keV) remains at a relatively 
constant flux until phase 0.9. Up to this point the flux at high 
energies appears to be flnely balanced between the intrinsic 
softening of the spectrum and the increasing luminosity as 
the stars approach each other. However, this balancing act is 
over by phase 0.9, after which the spectral softening rapidly 
accelerates. Between phase 0.9 and 1.0 (periastron) there is a 
precipitous collapse in the hard X-ray emission, as the stars 
move deep within the acceleration zone of the other's wind. 
The intrinsic emission is of comparable softness at phase 1.0 
and 1.1 (but brighter at phase 1.0 due to the reduced stel- 
lar separation), and is markedly harder by phase 0.2 (see 
Fig. 117b ') . These variations in the intrinsic spectra explain 
the intrinsic lightcurves shown in the top panels of Fig. [15] 

Fig. [TS] reveals that there is a strong hysteresis to the 
intrinsic emission, with large differences in the spectrum at 
identical stellar separations depending on whether the stars 
are approaching or receding from each other. The emission is 
much harder as the stars move together compared to when 
they separate: this is a natural consequence of the higher 
pre-shock wind speeds that are attained prior to reductions 
in the stellar separation. As the stars approach each other 
the conditions in the WCR reflect, to some extent, the hot 
and rarefled plasma created at earlier orbital phases. Simi- 
larly, as the stars recede the downstream conditions in the 
WCR reflect still the lower preshock velocities at earlier or- 
bital phases, and in extremum the cold, dense plasma cre- 
ated during periastron passage. The hysteresis is largest 
nearest periastron, when changes in the pre-shock condi- 
tions are at their most rapid, and smallest near apastron 
when the rate of change in the stellar separation is most 
sedate. The observed hysteresis is also partly due to varia- 
tions in the relative wind speeds towards each star - when 
the stellar separation is decreasing, the stars (and thus also 
their winds) have a component of their orbital velocity di- 
rected towards each other, which augments the wind speeds 
in the centre of mass frame. The opposite effect occurs when 
the stars recede from each other. This mechanism enhances 
the post-shock temperature in the WCR after apastron, and 
reduces it at comparable orbital phases prior to apastron. 
Hysteresis of the thermal radio and sub-mm emission also 
occurs (see Paper II). 

The hysteresis of the intrinsic emission is also clear when 
the luminosities are plotted against stellar separation, as 
shown in Fig. [I9ja)-(c). Interestingly, the intrinsic emission 
in the 2.5 — lOkeV band is stronger as the stars approach 
periastron, whereas the emission in the 0.1 — 0.5 keV band 
is stronger as the stars recede. The hard emission requires 
high temperature plasma, which the WCR is full of in the 
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Figure 15. X-ray lightcurves for model cwb4. The lightcurves were calculated over the energy bands 0.1 — 0.5 keV (left), 0.5 — 2.5 keV 
(middle), and 2.5 — 10 keV (right). The top panels show the intrinsic luminosity (red curves), while the bottom panels show the 
attenuated (observed) luminosity for viewing angles in the orbital plane (i = 90°) with </) = 0°, 90°, and 180° (red, green and blue 
curves, respectively), and directly above/below the orbital plane (i = 0° , pink curves). In all panels curves displaying a 1/dscp variation 
(with arbitrary normalization) are over-plotted. 
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Figure 16. Intrinsic and attenuated X-ray spectra from model cwb4 for viewing angles i = 90° and = 0° at a) periastron, and b) 
apastron. 



second half of the orbit, but which is comparatively lacking 
in the first half of the orbit. In contrast, emission in the 
soft band does not require such high temperatures. Instead 
it is strongest when the postshock densities are high, such 
as during the formation of cold clumps within the WCR. 
Some of this emission will be from intermediate temperature 
interface regions where hot plasma surrounds cooler clumps. 
In comparison the intrinsic emission in the 0.5 — 2.5 keV band 
displays a transitional state: the emission is brighter in some 
parts of the orbit when the stars are receding, and fainter 
in other parts. 

Armed with an understanding of how the intrinsic X- 
ray emission varies with orbital phase, we now examine the 
attenuated emission. The attenuated lightcurves shown in 
the bottom panels of Fig. [15] display behaviour which de- 
pends on the viewing angle of the observer. The lightcurves 
for observers at i = 90° and = or 180° are almost 
identical, particularly in the hard 2.5 — lOkeV band. Such 
symmetry is expected, given the identical stellar parameters. 
The 2.5 — lOkeV attenuated lightcurves show the closest be- 
haviour to their intrinsic counterpart, highlighting the abil- 
ity of hard X-rays to stream through the circumstellar envi- 



ronment relatively unaffected by absorption. For this same 
reason the 2.5 — 10 keV attenuated lightcurves also show very 
little change with the viewing angle of the observer, with the 
largest difference occuring at phases 0.8 — 0.9 when there is 
more attenuation for an observer at i = 90° and <j> = 90° 
than for other orientations, because the stars are eclipsing 
the apex of the WCR at this time (see Paper I). 

The attenuated 0.5 — 2.5 keV lightcurves all peak at 
phase 0.95, irrespective of the orientation of the observer, 
in agreement with the timing of the peak in the intrinsic 
lightcurve. However, the height and shape of the maximum 
in the attenuated lightcurves is dependent on the orienta- 
tion. The greatest luminosity in the attenuated 0.5 — 2.5keV 
lightcurves occurs for an observer viewing the system face 
on (j = 0°). In contrast, the timing of the maximum in the 
attenuated 0.1 — 0.5 keV lightcurves is strongly dependent on 
the viewer's orientation, ranging from phase 0.9 for observers 
in the orbital plane at viewing angles oi (j) = Q° and 180°, 
to phase 0.0 (periastron) for a viewing angle oi (jy = 90°. 
For an observer viewing the system face on the maximum 
occurs at phase 0.95 - incidentally, this is also the highest 
maximum seen in the 0.1 — 0.5 keV lightcurves. These differ- 
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Figure 18. Intrinsic X-ray spectra from model cwb4 as a function of orbital phase. In each panel the stars are at identical separations, 
and either receding (phase < 0.5) or approaching (phase > 0.5). The hysteresis of the spectra is clearly evident. 



ences in the timing of the maxima reflect the propensity for 
soft X-rays to be attenuated by the circumstellar environ- 
ment, and the dependence of the strength of this attenuation 
on the orientation of the observer. The strong circumstellar 
absorption near periastron arises from the enhanced wind 
densities around the WCR due to the reduced stellar sepa- 
ration and pre-shock wind speeds. 

Broad minima which are roughly centered on apas- 
tron occur in most of the attenuated 0.1 — 0.5 keV and 
0.5 — 2.5 keV lightcurves, though there is a slight maximum 
at apastron in the 0.1 — 0.5 keV lightcurve for an observer 
with i = 90° and 4> — 90° (since lines-of-sight from the apex 
of the WCR initially pass through the low opacity WCR). 



The minima following periastron are generally deepest at 
phase 0.1. The shape of the minimum is also generally quite 
smooth, though the 0.1 — 0.5 keV lightcurves for i — 90° 
and = 0° and 180° are noticeable for showing more struc- 
ture (Fig. I15r b) shows the luminosity levelling out between 
phases 0.95 — 1.0, before falling more steeply between phases 
0.0-0.05). 

Panels (d)-(f) of Fig. ll9l show the attenuated luminosi- 
ties in the three energy bands for an observer at i = 90° 
as a function of orbital separation. The 2.5 — lOkeV emis- 
sion is most similar to its intrinsic counterpart (Fig. 119b ). 
again illustrating the relative ease at which the hard X-rays 
travel through the circumstellar environment. In contrast. 
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Figure 19. Intrinsic (top row) and attenuated (bottom row) X-ray luminosities from model cwb4 as a function of orbital separation. 
The energy bands 0.1 — 0.5 keV (left), 0.5 — 2.5 keV (middle), and 2.5 — lOkeV (right), are shown. The attenuated lightcurves arc for an 
observer located directly above the orbital plane (j = 0°). Some orbital phases are marked on the plots. 



the emission in the 0.1 — 0.5 keV band suffers severe attenua- 
tion, and this has a large impact on the shape of its hysteresis 
curve (compare Figs. [T9k and d). As we have already seen, 
the attenuation is particularly severe after periastron, when 
there is not much hot, low opacity, plasma in the WCR. 

The eccentric orbit means that the emission from model 
cwb4 at various times resembles that from models cwbl and 
cwb2. Careful examination reveals that the periastron spec- 
trum (Fig. 116b .) is slightly harder than the phase 0.0 spec- 
trum from model cwbl (cf. Fig. |4^), with the Fe K emis- 
sion visible in the former plot. This reflects the fact that 
the downstream flow in model cwb4 contains hotter plasma 
(which was shocked when the winds previously collided at a 
higher speed) than in model cwbl. That they are otherwise 
so similar reflects the fact that emission at the apex of the 
WCR (which responds much quicker to changing pre-shock 
conditions than emission from far downstream) dominates 
the total emission in this model as the post-shock gas at the 
apex of the WCR rapidly becomes extremely radiative. 

Likewise, there is a high degree of similarity between 
the apastron spectrum (Fig. 116b ) and the phase 0.5 spec- 
trum from model cwb2 (cf. Fig.|4};), the former being slightly 
softer. This is again consistent with the recent history of the 
WCR. Their likeness reflects the fact that in model cwb4, 
the rate of change in the steUar separation is at its most 
sedate at apastron. The dynamical timescale for flow out of 
the system is then short compared to the timescale for sig- 
nificant orbital change. At phase 0.5, tflow ^ v''°]2 ~ ^-"^ '^^ 
while the time for the orbital separation to change 10 per 
cent from its apastron value is 0.19 Porb = 1-2 d. 

Fig.[2nishows the effective circumstellar column and op- 
tical depth as a function of energy and orbital phase for an 
observer with i = 90° and = 0°. The columns at phase 0.0 
and 0.5 bear some similarity to those obtained from models 
cwbl and cwb2, in the same way that the attenuated spec- 
tra do. Thus, to first order there is rough agreement between 



the emission and absorption characteristics of an eccentric 
system at periastron and apastron and circular systems of 
identical stellar separation, at least for the region of param- 
eter space covered in these models. 

Fig. I20f a) also shows that there are extremely large 
phase-dependent variations in the energy-dependent col- 
umn. The variation of the column to the high energy (e.g. 
5keV) emission is largely due to changes in the degree of 
occultation to this emission, with high occultation at con- 
junction (phase 0.0 and 0.5), and lesser occultation near 
quadrature (phase 0.14 and 0.86). There is a severe decline 
in the column to the high energy emission (due largely to 
changes in the degree of occultation) between periastron and 
phase 0.1. At phase 0.2 the observer views the WCR apex 
through hot plasma further downstream (see Fig. 10 in Pa- 
per I), and the circumstellar column declines at all ener- 
gies. By phase 0.3 one of the stars is already positioning 
itself in front of parts of the apex of the WCR, and the 
column to the high energy emission increases from its value 
at phase 0.2. The high energy column further increases to 
a maximum near phase 0.5. The column to the high en- 
ergy emission eases again by phase 0.8, while at phase 0.9 
the observer again views the WCR apex through hot plasma 
further downstream, which results in the lowest effective col- 
umn and optical depth at all energies and phases. 

In contrast to the 4 orders of magnitude variation in the 
column at 10 keV, the column to the low energy emission is 
surprisingly steady during the majority of the orbit, which 
reflects the large volume from which this emission arises. 
However, there is again a significant reduction in the low 
energy circumstellar column at phase 0.9 due to the reasons 
previously given. 

The phase dependent variation in the optical depth 
shown in Fig. I20f b) to a large part reflects the changes in 
the circumstellar columns commented on above. The varia- 
tion in the optical depth as a function of phase spans the 
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Figure 20. a) "Effective" circumstellar column (a) and optical depth (b) from model cwb4 for an observer in the orbital plane (i = 90°) 
as a function of energy and orbital phase. Note that occultation can severely affect the value of the "effective" column. It is therefore 
not always an accurate reflection of the column that the observed X-rays pass through (see Sec. l3.1."3l for further details). 



range 0.8 - 7 at 0.1 keV, 0.1 - 1 at 1 keV, and < 10^^ - 0.2 
at lOkeV. 



3.4.3 Spectral fits 

Tableland Fig.[2l]show the results of spectral fits to "fake" 
Suzaku spectra with a nominal exposure time of 20 ksec gen- 
erated from model cwb4. The spectra in Fig. l21l were specif- 
ically chosen to highlight the large spectral variations which 
occur over the course of the stellar orbit. The spectrum at 
phase 0.0 is almost as soft as the emission gets in this model 
(the spectrum at phase 0.05 is marginally softer), reflecting 
the strong coohng of the plasma in the WCR at this phase. 
By phase 0.2 (not shown) the spectrum is noticeably harder 
(and more luminous). The spectrum at phase 0.6 is about as 
hard as the emission gets, and shows a prominent FeK line. 
At phase 0.95 the spectrum is at its most luminous, and is 
again softer, reflecting the lower pre-shock wind speeds at 
this phase. There are now not enough counts at high ener- 
gies to detect the Fe K line in the binned spectra, though it 
is of course seen in our theoretical spectra. 

The spectral variability shown in Fig. [21] is reflected in 
changes in the values of the flt parameters (see Table 
which are plotted in Fig. (22] The temperature of the hottest 
mekal component shows significant variation, changing from 
0.66tJ!;j3 keV at phase 0.05, to 2.53155:?? keV at phase 0.7. 
In addition, significant enhancements in the normalization of 
the components occur as periastron is approached. Between 
apastron and phase 0.95 the normalization of the hot compo- 
nent increases by a factor of 4.6, far above the corresponding 
1/daep value. Panel d) shows the combined luminosity of the 
three mekal components. A comparison with Fig. 1 151 reveals 
that the fits do a good job of recovering the phase variation 
of the observed and also the intrinsic luminosity. It is clear 
that the fits infer the presence of significant circumstellar 
absorption between phases 0.95 — 1.05, which is responsible 
for the large difference in the intrinsic and ISM corrected 
luminosities. 

It is also noteable that the normalization of the hot 
component dominates those of the cooler components from 
phase 0.4 to 0.95, while the normalization to the warm (i.e. 
the second) component dominates at phase 0.0 and 0.05. 
There is no need for additional circumstellar absorption to 



Table 4. Comparison of the intrinsic 0.5 — lOkeV luminosity of 
the theoretical spectra with that obtained from the spectral fits. 
All values are in units of 10^^ ergs~^ . The percentage ratio of the 
inferred to the actual intrinsic luminosity is shown in brackets for 
each case. 





Model cwbl 


Model cwb2 


Intrinsic luminosity 


3.11 


1.81 


Chandra phase 0.0 


1.69 (54%) 


2.10 (116%) 


Chandra phase 0.25 


2.71 (87%) 


1.90 (105%) 


Suzaku phase 0.0 


1.65 (53%) 


1.75 (97%) 


Suzaku phase 0.25 


2.59 (83%) 


1.90 (105%) 



the mekal components at phase 0.9, which is consistent with 
the low value of the effective circumstellar column at this 
phase (see Fig. 120b .). While the fits do require substan- 
tial additional column to the warm component at phase 
0.0 (1.13lo;2i ^ lO^^cm"^), this extra absorption is about 
3 times higher than the effective circumstellar column at 
0.49 keV, as shown in Fig. 120b .). It is also puzzling why the 
spectral fitting did not require extra absorption (above the 
ISM value) to the cold and hot mekal components at this 
phase, despite these making significant contributions to the 
overall observed emission (see the top left panel in Fig. I2ip . 



4 COMPARISON TO OTHER NUMERICAL 
MODELS 

The X-ray emission from 0-|-0-star CWBs has 
been investigated using fully hydrodyna mical models 
(|Pittard fc Stovons"l997': 'Pittard et al]|2000 l) and "hydrid" 
models dAntokhin ct al. 2004; Parkin fc Pittard 20081). The 



analysis in the current work is a major improvement from 
that in iPittard fc Stevens! (|l997h . where the X-ray calcu- 
lations were based on 2D axisymmetric hydrodynamical 
models in which the winds instantaneously accelerated to 
their terminal speeds. As such, the plasma temperatures 
returned fr om model cwbl are much lower than those from 
model A in lPittard fc Stevens! (| 19971). The hydr o dynam ical 
model underlying the analysis in Pittard et al.l (|2000! ) did 
consider the radiatively-driven acceleration of the winds. 
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Figure 21. Three-temperature fits to "fake" Suzaku spectra from model cwb4 for an observer in the orbital plane viewing the system 
at conjunction [ip = 0°). a) Phase 0.0, b) phase 0.6, c) phase 0.95. The dramatic variation of the spectra with orbital phase is clearly 
shown. The assumed exposure time is 20 ksec in each case. 
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Figure 22. Orbital phase variability of the fit parameters to the "fake" Suzaku spectra from model cwb4 for an observer with i = 90° 
and (/) = 0°. Panel a) shows the temperature variation, panel b) the variation in the absorbing column (the displayed values include the 
assumed ISM column of A'^h = 10^^ cm^^), and panel c) the variation of the normalization, of each fit component. Panel d) shows the 
observed, ISM corrected and intrinsic luminosity of each component combined. 



but remained 2D (|Pittardl Il998l ). This has subsequently 
been shown to introduce an incorre ct phase dependence to 
the volume of hot eas in the WCR (|Lemaster et al.ll2007l '). 

In order to spatial ly resolve the post-sho ck cooling in 
highly radiative WCRs. P^ntokhin et al.l (|2004l ) introduced a 
model for the X-ray emission in which the local post-shock 
cooling was decoupled from a global solution of the ram- 
pressure balance. This work has the advantage of a low com- 
putational cost, allowing the rapid exploration of parameter 
space. However, its drawbacks include the neglect of orbital 
motion and the use of simple "beta" velocity-laws to mimic 
the winds' acceleration. The X-ray lightcurve of an equal 
winds system (model A) shows two minima per orbit, with 
dual symmetry about both conjunction and quadrature (this 
symmetry is broken in our models due to the aberration of 
the WCR). Strong absorption occurs when the observer's 
line-of-sight is down one of the arms of the WCR. This sig- 
nature is also seen from our model cwbl, though is greatly 



reduced in strength, because of the downstream curvature 
of the WCR. 

Three-dimensional colliding winds simulations of 
longer-period O+O systems (with orbital periods of 1 month 
and 1 year ), and eccentricity of 0.3 , have recently been pre- 
sented by IParkin fc PittardI (|2008l '). This work solved the 
equations of ram-pressure balance assuming that the winds 
collide at their terminal velocities, and adopted an analyt- 
ical expression for the abberation angle, to construct a 3D 
surface of the head of the WCR. Further downstream, the 
gas in the WCR was assumed to behave ballistically, which 
with the orbital motion of the stars results in a downstream 
curvature of the WCR. Again this approach benefits from its 
computational speed. The X-ray emission from this model 
was then calculated by mapping onto the surface of the 
WCR the emission calculated from a 2D hydrodynamical 
model. In all cases the WCR was assumed to be adiabatic. 
Radiative transfer through the computational volume with 
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Table 3. Spectral fitting results to simulated spectra from models cwbl, cwb2, and cwb3. The normalization of the mekal components 
(norm) is defined as [10~^^/(47rZ)2)] J ncnudV, where D is the distance to the source (in cm), and Ue and njj are the electron and 
hydrogen number densities (in cm"^)^ respectively. The indicated range in each fit parameter is the 90 per cent confidence interval for 
one interesting parameter. The last three columns give the observed fiux, the fiux corrected for ISM absorption, and the intrinsic flux, 
between 0.5 and 10 keV. All fits are to "fake" observations of 10 ksec duration, except the Suzaku observations of models cwb2 and cwb3, 
where the assumed exposure times are 40 ksec and 20 ksec, respectively. 



Model Phase Model fit 


kT 


Nn 


Norm 


Xl (d.o.f.) 


Observed fiux ISM corrected fiux Intrinsic fiux 




(kcV) 


(1022 cm-2) 


(10-3) 




(10-12 erg cm-2 g-i) 



cwbl 


0.0 


2T 


0.28g:ii 


lO-is 

'-'•^0-1 
590-67 


2 023-45 

4 074-58 
^■'^'2.86 


cwbl 


0.25 


2T 


250-27 

670-70 
'-'■"'0-64 


10-13 


4 iq5.27 

R r)Q6-32 
"-^^5-62 


cwb2 


0.0 


3T 


1 Olios 

2 473-25 


lO-is 

0.250-34 

0.3lg-f 


2.693;6| 
3.45|:f| 


cwb2 


0.25 


3T 


320-34 
840-93 


O.lOg-32 
10-24 

10-15 

'-'■^0-1 


9 1S5-67 

1 fiQ2.45 
l-DOl.44 
c c:tr5.84 
■^■"^'^5.09 



Chandra fits 

1.27 (93) 2.27 3.70 37.0 

2.21 2.74 10.4 

1.71 (110) 4.21 71.3 71.3 

5.31 6.81 15.5 

1.13 (165) 3.51 5.46 5.46 

4.01 4.90 6.99 

3.68 4.00 5.14 

1.04 (171) 2.75 4.31 4.34 

2.86 3.75 3.75 

6.65 7.76 7.76 

Suzaku fits 

cwbl 0.0 3T O.llg-J'* 0.17^-f 4.2l69g6 1.50 (324) 0.30 0.64 1.08 

0.3lg-g 0.15g-23 2.173^63 2.19 3.40 4.24 

0.69g-^| 0.57j]-62 3.343^65 gg 2.45 8.46 

cwbl 0.25 3T 0.120-14 0.IO2O-21 3.17^-^1 2.28 (401) 0.58 1.26 1.28 

0.3l[j-|J O.I4O-14 4.645^47 4 85 7 57 g 03 

0.740-74 0.410-43 4.604-|0 3.91 4.89 11.4 

cwb2 0.0 3T 0.33g-|f O-'^o i* l-S^i'tg ^-^^ (295) 2.46 3.82 3.82 

0.88g-|| O.I7O-28 2.23f-?| 3.26 4.16 4.96 

2.38|-f^ O-Iq?^ 4.10|-|^ 5.08 5.84 5.84 

cwb2 0.25 3T 0.3lg-g O.lfj'P 2.172|i 1.17 (319) 2.70 4.26 4.26 

0.83g-^^ O.I3O-31 1.9lf-J| 3.15 4.10 4.41 

2.142^26 0.1{;-17 5.16|-|g 6.24 7.25 7.25 

cwb3 0.0 3T 0.30g-|^ O-Iq P l-82f-gf 1.21 (339) 2.20 3.50 3.50 

0.84[j-|^ 0.15ji-f-* l-7lfi3 2.65 3.43 3.93 

2.092-21 0.lj;-f 2.122|6 2.55 2.97 2.97 

cwb3 0.25 3T 0.30S-|i 0.1°-^ 1-69H| 1.35 (353) 2.02 3.23 3.23 

0.85g-|f O.iof^l'^ I-SOJIt 2.62 3.43 3.23 

1.79j-ff 0.1H2 2.60ill 3.08 3.63 3.63 

cwb3 0.5 3T 0.29[j-|| O.lJJ-p l-75j-|:j: 1.14 (324) 2.05 3.30 3.30 

0.850|| O.I4O-24 1.572-09 2.48 3.21 3.60 

1.66j|f 0.46j|-^f 2.67^-?J 2.14 2.35 3.84 

cwb3 0.75 3T 0.30g-|J O.ljj-p l-75?'g3 0.94 (356) 2.12 3.38 3.38 

0.85g|| O.ljj^* l-66j|| 2.89 3.78 3.78 

1.992-08 O.lg-12 2.5l|-|g 3.02 3.52 3.52 



the inclusion of opacity then allowed the generation of syn- 
thetic spectra and lightcurves. It was found that the X-ray 
flux generally followed a 1/dscp scaling, but could also dis- 
play absorption related variations. While in principle the 
inclination and orientation of the system can be constrained 
from the shape of the X-ray lightcurve, in practice the lack 
of significant circumstellar absorption in O-f O star systems 
with periods of order 1 yr will make this very difficult. The 
situation of course improves in shorter period systems, and 
in systems with reasonably disparate wind densities (e.g. 
WR-FO, LBV-FO, and LBV-^WR systems). 



5 COMPARISON TO OBSERVATIONS 

In this section the results from our models are compared 
to recent X-ray observations of short period 0-|-0-star sys- 
tems. We focus first on systems where the WCR is expected 
to be highly radiative, then on systems where the WCR is 
expected to be more adiabatic, and then on systems with 
eccentric orbits and unequal winds. We finish by discussing 
the observational evidence for non-thermal X-ray emission. 



5.1 Highly radiative systems 

Model cwbl is similar to HD 215835 (DHCep), HD 165052, 
and HD 159176. The X-ray emission from each of these sys- 
tems is discussed below. 
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Table 5. As Tabic [3l but for model cwb4. All the fits arc of 3-tcmperature mekal models to simulated Suzaku spectra of 20ksec duration, 
for an assumed inclination i = 90° . 



Phase 


(keV) 


AT 

(10^'^ cm ^) 


Norm 

(10 3) 


^,2 / J ^ r \ 
(d.o.t.) 


Observed flux 


ISM corrected flux 
(10 ^-^ergcm ^s ^) 


Intrinsic flux 


0.00 




O.lO-Ji 




1.51 (253) 


6.04 


10.0 


10.0 




0.49{]i| 


1.23i-« 


1 q -118.4 
^•^-■'9.52 




1.72 


1.99 


33.7 




0.820-85 

^ " 0.79 


0.10-16 


1 qq2.25 
1.62 




3.33 


4.38 


4.38 


0.05 


0.262-31 


O.1O-" 

0.1 


Qi 3.53 

Z.Olo Act 


1.41 (218) 


2.94 


4.91 


4.91 




0.30 


1.24} 

1.04 


18.8ff-5 

11.4 




1.39 


1.64 


41.5 




0.66°-" 

Q-Qo 


0.10-20 

0.1 


1.862-59 




3.48 


4.76 


4.76 


0.10 


0.210-25 

U. LO 


0.10-^5 


3 75 
2.352 02 


1.32 (209) 


1.91 


3.48 


3.48 




O.47O-54 


0.10:1" 


l-78?il 




3.07 


4.44 


4.44 




O.8OH? 

0.75 


O.43O-5O 
0.37 


4 1/ 
Z.68 




2.90 


3.57 


8.30 


0.20 


0.290-28 




2.973-29 


1.35 (290) 


3.51 


5.62 


5.62 




0.800-f2 


0.33q 2g 


6.03|:^| 




6.23 


7.83 


14.3 




l-63}i8 

1.48 


O.lO-f 


1.752:11 

1.41 




2.13 


2.53 


2.53 


0.30 


0.30S-^J 

^^^0.28 


^0.1 


2.452-91 

2.23 


1.12 (298) 


2.94 


4.69 


4.69 




0.830-86 


0.370-43 

U.oU 


4.845;75 




4.59 


5.68 


11.2 




2.02f|° 

1.86 


0.10-20 

0.1 


3.032 gg 
2.10?-gt 




3.64 


4.24 


4.24 


0.40 


0.3ip| 


O.io-f 


1.05 (287) 


2.61 


4.12 


4.12 




0.88n'ec 

(J.QO 


0.240-35 


3.114.12 




3.91 


4.92 


6.96 




2.022;3| 


0.27o:f 


3 094:86 

O.U4 




3.86 


4.28 


5.49 


0.50 


. 3 T r\ 1-1 Q 


O.1O-19 


0.651-11 


1.02 (267) 


0.94 


1.42 


1.42 






0.180-f 


2.333:06 




3.33 


4.27 


5.32 




2.28?-39 

1.92 


O.1O-35 

0.1 


3.82t90 

3.52 




4.68 


5.40 


5.40 


0.60 




0.10-^5 


0.620-93 


1.02 (267) 


0.73 


1.17 


1.17 




o.87o:|o 


0.10:19 


1 702.04 
1.63 




3.09 


4.03 


4.03 






0.290-f 


c 045.59 




4.87 


5.40 


7.04 


0.70 


n 640-70 

^•°^0.38 


lO-is 


i-i7j:|| 


1.00 (303) 


2.19 


3.01 


3.01 




^•■^^0.86 


SI 0-48 


2.683-40 




2.77 


3.33 


5.34 




2 ^■^i.rs 

^•'-'"'2.36 


10-14 


4-231:25 




5.35 


6.12 


6.12 


0.80 


n S50.40 


l0-13 


1 431-^3 

l-^'^l.lO 


1.14 (378) 


1.99 


3.04 


3.04 




n S^O.89 
U.80(j g2 


270-34 


3-85t:i 




4.47 


5.60 


8.77 




2.082-15 


10-12 






8.63 


10.0 


10.0 


0.90 


0.33o;34 


lO-ii 


4.04|;45 


1.24 (490) 


5.30 


8.23 


8.23 






10-14 


c 006.66 
"^•^^5.64 




10.3 


13.3 


13.3 




1 741-78 

'J -^^0.31 


lO-ii 


n 912-3 
11-^11.0 




14.2 


16.8 


16.8 


0.95 


10-12 


6 146-94 
"•1^5.34 


1.47 (427) 


8.24 


12.7 


12.7 




0-738:6^? 


10-12 


"•^^5.68 




12.2 


16.4 


16.4 




0-981:00 


610-72 


17 618-9 
l'-°16.1 




11.5 


13.3 


36.5 



5.1.1 HD 165052 

HD 165052 is a double-lined bi nary which has re cently been 
classified as 06.5V + 07.5V (| Arias et al.ll2002l ). The orbit 
is slightly eccentric (e=0.09), with an inclination thought to 
be arou nd i = 20° . This system was first considered as a 
CWB bv lLuo et al] (|l990h . The nature of the WCR remains 
unclear. An analysis of the wind ram pressures reveals that 
there may be no balance at periastron, though this might be 
possible at apastron. The shocked winds are likely to both 
be radiative (x ~ !)• 

A ROSAT lightcurve of HD 165052 presented by 
ICorcoranI (1 19961) was repha sed to the correct 2.95 d orbital 
period bv lArias et al.l (|2002h . The lightcurve shows two min- 
ima per orbit, roughly centered on phases 0.0 and 0.5 (con- 
junction), in good agreement with the lightcurves presented 
in Fig. [3]from model cwbl. However, the amplitude of vari- 
ation of the X-ray flux is a factor of 2 or so, which is slightly 
greater than the variation in the i = 30° lightcurves of 
model cwbl in Fig. |3] This difference may be caused by the 



slight orbital eccentricity, and/or may reflect the absence of 
a wind-wind balance or the unequal wind strengths. 

Unfortunately, HD 165052 lies just outside the field-of- 
view of recent X-ray observations: a serendipitous XMM- 
Newton observation of the Lago on Nebula (M8) an d the very 
young open cluster NGC6530 (|Rauw et al.ll2002c) centered 
on the non-thermal radio emitting O-star 9 Sgr Jllauw et al.l 
2002bl ). and Chandra observations of NGC6530 ( Damiariil 
200J) and 9 Sgr (PI Gagne, not published). Dedicated ob- 



servations of HD 165052 are clearly necessary. 



5.1.2 HD 159176 

A single observation of the 07V-f07V system 
HD 15 9176 was made with XMM-Newton in March 
2001 (|De Becker et all l2004bl '). HD 159176 has an orbital 
period of 3.367 d, and a suspected inclination angle i ~ 50° . 
Since the stars (and winds) are identical the winds should 
collide at the system centre of mass, half-way between 
the stars. The collision should be radiative, even if the 
mass-loss rates are an order of magnitude lower than the 
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Figure 23. Two-tompcraturc mckal fit to combined XMM- 
Newton PN and MOSl+2 spectra of DH Cep. An interstellar col- 
umn of A^H = 0.667 X 10^^ cm~^ has been assumed. 

values quoted in IPachoulaki^ (119961 ). The XMM-Newton 
observation was taken just after quadrature, corresponding 
to phases 0.28 — 0.41 (or, alternatively, phase 0.78 — 0.91) in 
our models (note that the phase reported in De Becker et 
al. is 0.53 — 0.66, but in their paper phase 0.0 corresponds 
to the maximum radial velocity of the primary compo- 
nent). The observed luminosity in the 0.4 — lOkeV band 
is 1.5 X 10'^'' ergs~^. Fig. [3] shows that the attenuated 
luminosity near quadrature in model cwbl is, in fact, very 
similar to that observed from HD 159176. T his is important, 
because theoretical models presented by iDe Becker et al.l 
(|2004bl ) which decoupled the local postshock cooling 
from the global dynamics of the WCR overpredicted the 
luminosity by at least a factor of 4. 

A two-temperature mekal fit to the combined EPIC 
and RGS data yielded kTi = 0.21 ± 0.02 keV and kT2 = 
0.96 ± 0.01 keV, with a circumstellar column Nh = 0.40 ± 
0.01 X 10^^ cm~^. The latter is in good agreement with the 
energy dependent column at quadrature from model cwbl 
as shown in Fig. Efa) . We further find that for an inclination 
angle i = 90°, at phase 0.25 the characteristic mekal temper- 
atures from spectral fits to model cwbl are 0.25lo keV and 
0.67±0.03 keV, and the circumstellar column to the hot com- 
ponent is 0.29+° g4 X 10 cm . Therefore, the agreement 
is quite good. The slightly higher temperature of the hot 
component in the fit to the real data may indicate slightly 
higher pre-shock wind speeds than obtained in our model. 

5.1.3 HD 215835 (DH Cep ) 

HD 215835 (DHCep), located in the young open clus- 
ter NGC7380, is an extremely tight 06V-I-07V bi- 
nary with a circular orbit of just 2.1109 d period 
I Penny. Gies fc Bagnuolclll997l ). The stars are tidally dis- 
torted, but at present are fully detached. The inclination 
remains slightly uncertain (35° < i < 51°). The significant 
optical polarization is consistent with scattering within a 
WCR, while its short-te rm variability m ay be related to in- 
stabilities in the WCR (|Corcoranlll99lf ). An analysis of the 
wind ram pressures reveals that there is unlikely to be any 
balance, and the primary wind should directly impact the 
secondary star. However, once again a detailed simulation 



of this system will be necessary to place this hypothesis on 
a firmer footing. 

DH Cep was observed by XMM-Newton on 
ig**" December, 2003, for « 30ksec with the MOSl 
and M0S2 cameras, and f» 28ksec with the PN camera. 
It was reduced using standard scripts and version 9.0.0 of 
SAS, and the spectrum was rebinned to obtain a minimum 
of 25 counts per bin. Fig. [23]shows a two-temperature mekal 
fit, with an assumed ISM column of 0.667 x 10^^ cm"^, the 
weighted average returned by the NH tool at the HEASARC 
website. The fit is good (xl ~ 1-18), and has the following 
parameters: kTi = 0.24l° °J keV, kT2 = 0.761551^^ keV, 
iVHi = 0.671°;°^ X 10^2 gm-^ = ^-O^toZ x 10^^cm-^ 
normi = 2.07+o;24 x 10"^ and norm2 = 4,.29to.fl x 10"^. 
The returned temperatures are similar to those from 
model cwbl. The fit needed no extra circumstellar ab- 
sorption to the soft component (as also found for model 
cwbl), but favoured a circumstellar column of about 
0.4 X 10'^'^ cm~'^ to the hard component (again in good 
agreemen t with model cwbl). Assuming a distan ce of 
3.73 kpc (|Massev. Johnson fc Degioia-EastwoodI Il995^ . the 
ISM corrected 0.5 - 10 keV luminosity i s 6.3 x 10^^ ergs"^ . 
With logLboi = 5.85 (|Pennv et al.1 Il997f l. we obtain 
ix/iboi = 2.3 X 10~^. The X-ray luminosity and Lx/Lboi 
values are slightly higher than obtained from model cwbl. 
It is entirely plausible that such a high X-ray luminosity 
could be obtained even if the primary wind crushes that 
of its companion, because in such a scenario the secondary 
star would nonetheless intercept a large fraction of the 
primary's wind due to its proximity. 

The slightly hotter temperature of the 2nd mekal 
component compared to model cwbl may indicate that 
the primary wind indeed collides directly with the sur- 
face of the companion star, since it then has more room 
to accelerate before its collision, and thus the ability 
(ignoring potential braking and inhibi tion effects - see 
iGavlev. Owocki fc Cranmeii (|l997l ) and IStevens fc PoUockl 
(|1994| ). respectively) to reach higher pre-shock speeds. 

5.2 Intermediate and adiabatic systems 

Systems with slightly longer orbital periods and wider stellar 
separations, some of which bear similarities to model cwb2, 
are now considered. Specifically, observations of HD 93161A 
and HD 47129 (Plaskett's star) are discussed. 

5.2.1 HD 931 61 A 

HD 93161A is an 08V + 09V system with an orbital period 
of 8.566 d. The mimimum masses from the orbital solution 
are quite large, which suggests that the inclination of the 
orbit is high, most probably > 75° (|Naze et al.1 120051 ). The 
winds achieve a ram pressure balance, and both shocked 
winds are largely adiabatic (the primary's more so: we esti- 
mate that xi ~ 30 and X2 ~ 5) . 

The system was observed five times with XMM-Newton 
during 2000 and 2001, but no significant variations at ener- 
gies above 0.5 keV were detected. This is nicely consistent 
with the theoretical lightcurves from model cwb2 shown in 
Fig. |3l which display little orbital variation in energy bands 
above 0.5 keV due to; i) the low attenuation through the 
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circumstellar winds caused by the relatively low mass-loss 
rates and the relatively wide stellar separation, and ii) low 
occultation due to the large size of the WCR relative to the 
stars. 

The observed luminosity in the 0.4 — lOkeV band is 
1.2 X lO^'^ ergs~^ , which is an order of magnitude fainter 
than model cwb2. At first glance this would appear to be 
a big problem, but there are two important factors which 
can reconcile this issue. Firstly, the interstellar column to 
this system (A^^h = 4.5 x 10^^ cm^^) is significantly higher 
than assumed for our models. Secondly, the winds in this 
system are likely to be much more feeble compared to the 
06V winds assumed in our model. Since the X-ray luminos- 
ity Lx oc M in adiabatic systems, a reduction of a factor of 
3 or so in the mass-loss rates would bring the observed and 
theoretical luminosities into better agreement. Such a reduc- 
tion is consistent with the expected change in the mass-loss 
rate between an 06 a nd an 08/9 main-sequence star. 

iNaze et all (|2005l ) find that two-temperature mekal fits 
to the X-ray spectra of HD 93161A yield average values of 
fcTi = 0.28 ± 0.02 keV and kT2 = 0.76 ± 0.17 keV. The de- 
rived circumstellar column varies between 0.45 x 10^^ cm~^ 
(at phases 0.18 - 0.5 3) and 0.86 x 10^^ cm ~^ (at phase 
0.75). More recently, lAntokhin et al.l (|2008l ) performed a 
2-temperature fit to the combined data, finding tempera- 
tures of 0.27t|5;E!3 keV and 1.16l[5:g keV. The hotter tem- 
perature from this fit is greater than the temperature of 
the hot c omponent in any o f the 5 individual spectra an- 
alyzed by iNaze et all (|2005l ). This highUghts some of the 
non- uniqueness issues of spectral fits to medium-resolution 
spectra which is well documented in the literature. 

In comparison, two-temperature mekal fits to model 
cwb2 are generally pretty poor (typically xt > 2), with 
three-temperatures needed to provide satisfactory fits. With 
this in mind, typical temperatures obtained from two- 
component fits are 0.25 keV and 1.4 keV, with the analysis 
of simulated Chandra spectra returning higher temperatures 
for the hotter component (1.64 keV and 1.68 keV, at phase 
0.0 and 0.25, respectively) than the corresponding analysis 
of a "fake" Suzaku dataset (1.07 keV and 1.35 keV, at phase 
0.0 and 0.25, respectively). Excluding one of the Suzaku 
analyzes, the values of kT2 considerably exceed those re- 
ported from HD 9316 1A, even the higher value found by 
lAntokhin et al.l (|2008l ). However, this is not too surprising 
given that the wind speeds of the 08V and 09V stars in 
HD 93161 A are likely to be considerably lower than those 
assumed from the 06V stars in model cwb2. 



5.2.2 HD 47129 (Plaskett's star) 

X-ray observations of systems containing evolved O-stars 
have also been presented in the recent literature. One 
notable analysis concerns Plaskett's star (HD 47129), an 
08HI/I + 07.5HI system wi th a circular orbit of pe- 
riod 14.4 d (jLinder et all l2008h . The stars are very mas- 
sive (54 Mq and 56 Mq for the primary and secondary 
star respecti v ely, a ssuming an orbital inclination of 71°). 
iLinder et al.l (|2008h note that the secondary star is de- 
formed due to its large rotational velocity, which results 
in a non-uniform temperature distribution, and speculate 
that the wind is confined near the equatorial plane. In this 



way the apex of the WCR can occur closer to the pri- 
mary star (even if its wind has the overall greater mass-loss 
rate), as suggested from the analysis of o ptical emission lines 
IWiggs fc Gieslll992l : iLinder et al.]|2008l ). 

XMM -Newton obse r vation s of Plaskett's star were pre- 
sented by ILinder et all (|2006D . While there appears to be 
a minimum in the observed count rate when the primary 
star is in front, the amplitude of variability remains uncer- 
tain, as it is quite small between two XMM-Newton obser- 
vations, but larger in archival ROSAT HRI observations. 
There is a pressing need for additional data to resolve this 
issue, and to examine whether there is an additional mini- 
mum during eac h orbit. The interstellar absorbing column is 
1.5 X 10^^ cm"^ jDiplas fc Savage|[l993 ) . similar to the value 
assumed for our models. The observed 0.5 — 10 keV luminos- 
ity from the XMM-Newton observations (6.9 x 10^^ erg s~^ ) 
is also similar to the luminosities obtained in models cwb2 
and cwb3, while a minimum in the luminosity near phase 
0.0 is again consistent with the models. 

There is some evidence (e.g. a very faint FeK line) 
that the hot X-ray emitting plasma may not be in ther- 
mal equilibrium, or that there is a power-law tail to the 
hard emission. While the latter may indicate that non- 
thermal processes, such as particle acceleration, are oc- 
curing, it could also be an artifact of the fitting method: 
additional fits using a bremsstrahlung model to fit the 
continuum plus individual delta functions to fit the lines 
did not require a power-law component. The tempera- 
ture returned for the bremsstrahlung component is 2.22 ± 
0.1 keV. Three-component mekal fits return temperatures of 
0.3llEJ:^i keV, 0.74lg;°^ keV and 2.42to.?6 keV. No extra col- 
umn was needed to the softest component, but circumstellar 
columns of 4.7155;? x 10^^ cm~^ and 'i.OtYi x 10^^ cm~^ were 
needed for the intermediate and hard components, respec- 
tively. 

The temperatures reported above are comparable to 
those obtained by three-temperature mekal fits to the syn- 
thetic spectra generated from models cwb2 and cwb3, and 
noted in Table |3| In addition, a lower column to the harder 
than to the intermediate component is consistent with the 
phase 0.25 results in Fig. [3e). However, a specific model of 
this system is needed to determine whether all the observa- 
tional data can be accurately reproduced. 

5.3 Eccentric/unequal wind systems 

We now discuss systems with eccentric orbits and/or un- 
equal winds, bearing some similarities to models cwb3 and 
cwb4. 

5.3.1 HD 93403 

Four XMM-Newton observations of HD 93403, an 
05.51 -I-07V binary with a 15.093 d orbit of eccentric- 
ity e = 0.234, reveal flux variations of 20 per cent amplitude 
in the 0.5 — 2.5 keV band, with a m inimum centered o n 
apastron and a peak at periastron (|Rauw et al.l l2002al ). 
In the softer 0.5 — 1.0 keV band the minimum occurs at 
phase 0.75, when the denser wind of the 05.51 primary is 
in front. Va riations in the harde r 2.5 — lOkeV band are less 
significant. iRauw et al] (|2002al ) found that HD 93403 ap- 
pears less overluminous in X-rays than previously thought. 
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and suggested that a significant fraction of tlie total X-ray 
emission may arise in line-instability shocks intrinsic to the 
winds. The orbital inclination is thought to be ~ 30°. 

The spectrum obtained nearest to periastron is not 
markedly softer than the spectra obtained at other phases. A 
two-temperature mekal fit yields temperatures of ~ 0.23 keV 
and 1.06 keV, with the harder component showing more vari- 
ability. The observed luminosity (assuming a distance of 
3.2 kpc) in the 0.5 - lOkeV band is « 1.6 x 10^^ erg s"^ 
The circumstellar column was found to be ~ 8 x 10^" cm^^ 
and ~ 4.6 x lO^'^ cm~^ to the soft and hard mekal compo- 
nents, respectively. 

This system combines features of both models cwb3 and 
cwb4 (unequal winds and an eccentric orbit). However, un- 
like model cwb4, analysis of the wind ram pressures suggests 
that the winds may not be able to sustain a stable balance 
at periastron, with the result that the more powerful wind 
from the 05.5I star overwhelms that of its companion, and 
crashes directly onto its surface. Having said this, a full ex- 
amination of the potential effects of radiative inhibition and 
braking is necessary to strengthen this hypothesis. Assum- 
ing that the primary wind does impact the companion star 
at periastron, it will be shocked to high temperatures, and 
should radiate reasonably efficiently (x ^ 1)- A stable wind 
balance should be possible at apastron, but whether the 07V 
star can reestablish its wind towards its larger companion 
also remains to be seen. If a wind balance occurs, we esti- 
mate that both the shocked primary and secondary winds 
should be marginally adiabatic (xi ~ 2.5, X2 ~ 5). 

Bearing in mind that the nature of the WCR in this 
system could be quite different to that in models cwb3 and 
cwb4, we compare in the following the X-ray properties 
of this system with those from our models. We first note 
that the observed luminosity is approximately a factor of 2 
greater than the luminosity from model cwb3, but is compa- 
rable to the luminosity from model cwb4 over most of the or- 
bit, though the interstellar absorption column to HD 93403 
at 3.67 X 10^^ cm~^ is higher than that assumed in our mod- 
els. The circumstellar columns to the mekal components 
are comparable to the energy-dependent columns shown in 
Fig. [SJe). We can also make a direct comparison to the re- 
sults of a two-temperature mekal fit to a simulated lOksec 
Suzaku observation of model cwb3 at i = 90° and phase 
0.0. Though such a fit is poor {xl = 1-81, with 224 d.o.f.), 
we find that the fit temperatures (fcTi = 0.26 ± 0.02 keV, 
fcTb = 1.21 ± 0.05 keV) are nevertheless roughly compara- 
ble to the values obtained from analysis of the observations 
of HD 93403. The circumstellar column to the cooler com- 
ponent is not in very good agreement (A^Hi = 0.231^q f)? ^ 
lO^'^ cm~^), though the column to the hotter component is 
a better match (Ah^ = 0.31 ± 0.08 x 10^^ cm"^). 

5.3.2 Cyg OB2#8A 

Cyg OB2#8A is a longer period (21.9 d) 06If + 
05.5in(f) binary with an orbital e ccentricity e = 0.24 
(|De Becker. Rauw fc ManfroidI 120041 '). It has long been 
known as a very bright X-ray source. XMM-Newton obser- 
vations, together with archive RO SAT and ASGA data, re- 
veal that the observed emission reaches a maximum near 
phase 0.75, with a minimum likely to oc cur shortly after pe- 
riastron passage (|De Becker et al]|2006h . The intrinsic X-ray 



luminosity varies between 1.0 — 1.9 x 10'^* ergs~^ , giving 
an excess of ~ 20 compared to the canonical relationship for 
single star^. Clearly, this is caused by the collision of two 
very powerful winds. The orbital modulation of the emission 
matches roughly the 1/dsep variation expected if the winds 
are adiabatic. The emission is best fit by a 3-temperature 
mekal model with plasma temperatures of approximately 
0.24 keV, 0.80 keV, and 1.76 keV, though the hottest com- 
ponent shows a significant decline in temperature (by about 
15 per cen t) between apastron an d periastron. 

While iDe Becker et al.l (|2006l ') estimated that the post- 
shock winds would be radiative around the entire orbit, 
the assumed mass-loss rates may have been too high, 
and the terminal velocities too low. With mass-loss rates 
and terminal velocities more ty pical of recent determina- 
tions for the st ellar types (e.g. iRepolust. Puis fc Herrerol 
I2OO4I : lMartins~c t al. 2005a|) , we estimate that the primary 
wind possibly becomes radiative at periastron though is 
marginally adiabatic at apastron, while the secondary wind 
should remain adiabatic throughout the orbit. 

Since the nature of the WCR in Cyg OB2#8A is clearly 
different to that in any of the models calculated in Paper I, 
it is unsurprising that the behaviour and characteristics of 
its X-ray emission do not closely match those from any of 
our models. Nevertheless, it is interesting to compare this 
system to our model cwb4 which also displays variations in 
its fit parameters (see Tableland Fig. I22p. In model cwb4 
the hot component from 3-temperature fits declines by 65 
per cent, a much larger fall than seen from the actual ob- 
servations of Cyg OB2#8A. That the variability is much 
lower in Cyg OB2#8A is likely due to smaller changes in 
th e pre-shock wind spee ds around the orbit (see Table 11 
in lDe Becker et aLll2006l ) compared to model cwb4. In addi- 
tion, the minimum temperature of the hot component in the 
fits to model cwb4 actually occurs after periastron (phase 
0.05), this "time delay" being due to the history of the 
plasma temperature within the WCR. It would be interest- 
ing to obtain further X-ray observations of Cyg OB2:7^8A, 
with better phase coverage, to search for a similar effect, al- 
though the reduced eccentricity and longer orbital period 
may shorten or prevent a similar time delay from being 
seen. A specific hydrodynamical model of this system, which 
would allow further comparison of its emission in both the 
X-ray and radio domains (Cyg OB2#8A is also a strongly 
variable non-thermal radio source), would clearly be of in- 
terest. 



5.3.3 HD 37043 (i On) 

tOri is a highly eccentric (e = 0.764), 09III -I-BIIII b i- 
nary with a 29.134 d orbital period (|Marchenko et aLlbOOGD . 
which may have formed via an exch ange of binary compo - 
nents in a binary-binary collision (|Bagnuolo et al.l I2OOII ). 
It was observed twice by ASCA during 1997, with the 
obser vations timed to c oincide with periastron and apas- 
tron (jPittard et al.ll200d V The X-ray emission is bright be- 

^ Note, however, that this estimate may need to be substan- 
tially revised downwards if t he recent distance determination to 
Cyg OB2#5 of 925 ± 25 pc llLinder et al.ll2009l) is found to also 
be applicable to the other members of the association. 
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cause of its relative proximity [D ~ 450 pc). The observed 
0.5 — lOkeV luminosities are ~ 1.0 x 10'^^ ergs~^ in both 
pointings. Since the wind attenuation is negligible, and the 
interstellar column is low (A'h = 2 x lO^^cm"^), the in- 
trinsic X-ray luminosity of this system is also relatively low, 
consistent with the relatively feeble winds in this system. 
Surprisingly, there is no significant variation in the luminos- 
ity, spectral shape, or absorption between the periastron and 
apastron observations. Possible explanations for this lack of 
variability are that the WCR stays pinned to the surface of 
the secondary star throughout the entire orbit, or that in- 
trinsic shocks within the winds generated by the radiatively 
driven line de shadowing instability dominate the emission 
(|Pittardlll99i ). 

5.4 Non-thermal X-ray emission 

Since accelerated particles exist in many CWBs, the possibil- 
ity of inverse Compton X-ray emissi on has previou sly been 
raised by a number of authors (e.g. |Pollocklll987l ). Others 
have argued that non-thermal X-rays can be produced in the 
intrinsic wi nd shocks fo rmed through the line-deshadowing 
instability (|Chen &: Whi te 1991). A recent review of non- 
th ermal emis s ion pr ocesses in massive binaries can be found 
m iDe Beckeij (|2007l ). However non-thermal X-rays might be 
produced, their detection requires that their flux become sig- 
nificant relative to the thermal X-ray flux. Unfortunately, a 
campaign focused on O-stars with known non-thermal radio 
emission failed to unambiguously detect any non-thermal 
X-ray's at energ ies below lO keV (iRauw et al.|[2002bl . l2005l : 
iDe Becker et al.il2004d . |2005| . |2006| ). Onlv"from HD 159176 
was possible evidence o f a non-thermal component found 
(|De Becker et al.ll2004bh . As noted in Sec. 15. 1.21 this system 
has a very short period. The stellar winds likely collide at 
relatively low speeds, which limits the hardness of the ther- 
mal X-ray emission, allowing easier detection of a possible 
non-thermal component. 

The chance of unambiguously detecting non-thermal X- 
ray's increases with energy, as the emission from the thermal 
component falls off. It is partly for this reason that non- 
thermal X-ray emiss ion has been reporte d from a Suzaku 
observation of r] Car (|Sekie:uchi et al.ll2009^ . However, it ap- 
pears that, at least in the soft {E < lOkeV) domain, non- 
thermal X-rays are at best a relatively insignificant compo- 
nenlQ 



6 SUMMARY AND CONCLUSIONS 

This work investigates the X-ray emission from short-period 
0-|-0-star binaries where orbital effects and the acceleration 
of the winds are important. The structure and dynamics 
of the stellar winds and the wind-wind collision region for 
four distinct systems were previously calculated using a 3D 
hydrodynamical code (Paper I). We explore the emission 
arising from these models, under conditions where the WCR 
is radiative (model cwbl) or adiabatic (model cwb2), where 



* A summary of detections and upper limits o f CWB s at MeV, 
GcV and TeV energies can be found in lPittar 3 ll2009d) . 



the winds have unequal strengths (model cwb3), and where 
the orbit is eccentric (model cwb4). 

We find that model cwbl shows the greatest X-ray vari- 
ability of the three models with circular orbits (models cwbl, 
cwb2 and cwbS), due to the closer stellar separation and the 
higher densities surrounding the WCR. The lightcurves from 
each of these models are asymmetrical because of orbit in- 
duced aberration and curvature of the WCR. The asymme- 
try is greatest for model cwbl, which has the highest ratio 
of orbital to wind speeds, Uorb/fw, of the 3 simulations. The 
lightcurve asymmetry is always greatest for an observer lo- 
cated in the orbital plane {i = 90°), and generally decreases 
with decreasing i (as does the level of variability). Model 
cwb4 shows that the X-ray emission from systems with ec- 
centric orbits can be spectacularly variable, and can display 
a strong hysteresis around the orbit. In general, the emis- 
sion is softer while the stars are separating after periastron 
passage, and harder when the stars are approaching. The 
degree of hysteresis, and its duration, likely depends on the 
amplitude of changes to the ratio of the orbital period to 
the fiow time of the shocked plasma out of the system. As 
expected, the increase in the luminosity as periastron is ap- 
proached does not follow the 1/dsep scaling of an adiabatic, 
terminal wind speeds, WCR. The observed change also de- 
pends on the energy band and the observer's orientation. 
The ratio of the ISM corrected 0.5 — 10 keV X-ray luminos- 
ity to the system bolometric luminosity, -Lx/iboi, reaches 
values as high as 5.9 x lO^'' under favourable viewing an- 
gles, consistent with a very strong colliding winds signature. 
All-in-all, the diversity and richness of the simulated X-ray 
lightcurves matches well their observational counterparts. 

We use the difference between the intrinsic and attenu- 
ated spectra from the models to calculate "effective" cir- 
cumstellar absorbing columns which we demonstrate can 
be extremely energy dependent. Therefore, simple spectral 
fits which assume energy independent columns are over- 
simplified in the same way that modelling the spectrum 
with only a handful of different plasma temperatures is, and 
risk compromising the subsequent interpretation. We also 
fold our theoretical spectra through the response flies of the 
Chandra and Suzaku observatories, add poisson noise, and 
then analyze using the same method used for the analy- 
sis of real data. In many cases the flts return circumstel- 
lar columns in good agreement with the energy-dependent 
"effective" columns calculated directly from the theoretical 
spectra. However, there are also many instances where the 
returned columns have values which are not so close, and/or 
do not properly capture the trend of the effective circumstel- 
lar column with energy. We also flnd that despite our theo- 
retical spectra being calculated with solar abundances, ex- 
clusively sub-solar values are returned by the spectral fltting 
package to our simulated spectra when the global abundance 
is allowed to fit freely. The global abundances become closer 
to the true (solar) value as the complexity of the spectral 
model increases. These issues highlight some of the prob- 
lems associated with fitting overly simplified models to spec- 
tra, including non-uniqueness of the returned values, and of 
course are of wider signiflcance than just CWBs. Moreover, 
the intrinsic luminosities inferred from very short period sys- 
tems are often underestimated (by up to a factor of two), 
since the spectral models as used are unable to account for 
occultation losses. Despite these failings, in many cases the 
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resulting fit parameters are, at least to first order, similar to 
those obtained from the analysis of real observations. More 
stringent comparisons between theory and observations will 
require the modelling of specific systems, such as work cur- 
rently being conducted on 77 Car. 

In this paper we have examined the emission from 
main-sequence 0-|-0-star systems, in which circumstellar 
absorption is relatively minor. Circumstellar absorption be- 
comes more important as the stellar mass-loss rates in- 
crease, so that CWBs which contain a Wolf-Rayet (WR) 
star should display more significant orbital modulation of 
their lightcurves and spectra than seen in this work. For this 
reason, future work will explore the dynamics and emission 
from W R-I-Q systems, where dynamical effects like radiative 
braking (jCavlev et al.|[l997l ) will also be important. It will 
also be interesting to explore other parts of the CWB pa- 
rameter space. This could include: i) systems where a weaker 
wind is completely overwhelmed by a stronger wind (rele- 
vant to HD 93205, amongst others), ii) systems with an ec- 
centric orbit where all (or one side) of the shocked plasma 
in the WCR is radiative throughout the entire orbit (rele- 
vant to systems like HD 155248, an 07.5(f )III + Q7(f)III bi- 
nary with a 5.816 d orbit of eccentricity 0.127 (|Maver et al.l 
bOOll L which shows strong ph ase-locked X-ray emission 
with an asymmetric modulation (jSana et aLlbOO^ )). iii) sys- 
tems with more than two stars (relevant to HD 167971 and 
QZCar, amongst others). 

This work has also been limited to simulating 
medium resolution X-ray spectra. Future work will ex- 
amine X-ray line profiles, whi ch can directly probe the 
dynamics deep in the WCR (|Henlev. Stevens fc Pittaxdl 
12001 I2OO5I : iHenlev et ahl 120081 ). and intere sting physics 
such as non-equilibrium ionization (e.g. IPoUock et al.l 
I2OO5I) and non-equilibriu m electron and ion temperatures 
llZhekov fc Skinneij[2000f ). Further papers will examine the 
high energy non-thermal emission at X-ray and 7-ray ener- 
gies up to the TeV range. 

Despite this long to-do list, the present work has pro- 
vided further essential groundwork, highlighting the factors 
which affect various aspects of the X-ray emission from short 
period O-l-O star binaries. It lays the foundations for a bet- 
ter understanding of some of the key physics of such systems 
(e.g. radiatively driven winds, high Mach number shocks), 
which will be studied with future models of particular sys- 
tems. 
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